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A  Comparison  of  Moisture  Variables  in  the  Vertical 
Interpolations  of  a  4-D  Assimilation  System 


1.  INTRODUCTION 

Global  NWP  centers  routinely  generate  objective  analyses  of  the  global  atmos¬ 
phere  to  provide  an  accurate  depiction  of  its  structure  at  a  given  time  and  on  a 
uniform  three-dimensional  (3-D)  grid.  Characteristically,  the  analysis  function 
is  performed  in  the  context  of  a  global  four-dimensional  (4-D)  data  assimilation 
system.  The  crucial  components  of  a  4-D  assimilation  system  are  observations, 
objective  analysis  (interpolation  of  observations  to  a  uniform  grid),  initialization 
(balancing  of  mass  and  motion  fields),  and  short-range  forecast  (usually  of  6-  or 
12-hr  duration). 

NWP  centers  such  as  ECMWF,  *  NMC,^  and  AFGWC^  currently  use  similar 
4-D  data  assimilation  systems.*  All  three  centers  employ  6-hr  update  cycles,  as 
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Illustrated  in  Figure  1.  An  optimal  interpolation  (OI)  objective  analysis  of  the  ob¬ 
servations  updates  a  6 -hr  forecast  from  a  global  spectral  model  (GSM).  The  anal¬ 
yzed  mass  and  motion  fields  are  then  balanced  using  a  normal  mode  initialization 
(NMI)  procedure.  Next,  a  6-hr  GSM  forecast  provides  the  first-guess  fields  for 
the  next  OI  analysis.  Necessarily,  this  process  involves  vertical  interpolations  to 
transfer  the  analysis  fields  on  pressure  surfaces  to  the  o-coordinate  surfaces  of 
the  GSM  before  the  forecast  and  vice  versa  after  the  forecast.  These  interpolation 
steps  are  referred  to  as  the  preprocessor  and  postprocessor,  respectively. 
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Figure  1.  Schematic  Diagram  of  a  6-Hourly  Intermittent  Data  Assimilation 
System  (Top)  Contrasted  With  a  Standard  Forecast  Execution  Over  the  Same 
Time  Period  (Bottom) 

A  major  difficulty  of  data  assimilation  is  that  because  of  biases  and  inaccur¬ 
acies  in  the  forecast,  initialization,  and  pre- and  postprocessing,  continuous  rep¬ 
etitions  of  the  6-hr  cycle  contribute  an  undesirable  climatic  drift  to  the  assimila¬ 
tion.  This  climatic  drift  can  occur  despite  the  corrective  effects  of  the  objective 
analysis.  This  is  particularly  true  in  the  case  of  humidity  assimilation,  since 
characteristically  there  are  fewer  humidity  observations  than  wind  and  mass  ob¬ 
servations  to  support  the  objective  analysis.  Heretofore,  most  attention  to  the 
climate  drift  problem  has  focused  on  the  initialization  and  forecast  model  compo¬ 
nents,  with  lesser  attention  paid  to  the  vertical  interpolation  interfaces  in  pre- 
and  postprocessing.  Similarly,  even  less  attention  has  been  given  to  climatic 
drift  in  the  humidity  fields  as  compared  to  the  mass  and  motion  fields. 


Historically,  the  lack  of  attention  to  global  humidity  analysis  derives  from 
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studies  *  that  show  that  outside  the  tropics  large-scale  models  fairly  quickly 
generate  reasonable  hum  'ity  fields  after  one  or  two  days,  independent  of  the  in¬ 
put  humidity  analysis.  Rather,  the  humidity  fields  evolve  in  response  to  the 
large-scale  model's  developing  vertical  motion  fields.  Hence  in  the  extratropics, 
the  task  of  humidity  analysis  may  be  redundant  for  purposes  of  initializing  many 
large-scale  NWP  models.  In  this  context,  FNOC  chooses  not  to  perform  a  humid- 
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ity  analysis  in  its  global  data  assimilation  system.  Instead,  the  first-guess  fore¬ 
cast  field  of  humidity  becomes  the  initial  moisture  field  for  the  subsequent  assim¬ 
ilation  cycle. 

However,  objective  analysis  is  not  necessarily  performed  merely  to  initialize 
NWP  models.  Often,  objective  analysis  is  important  for  depicting  atmospheric 
structure  for  diagnostic  general  circulation  studies  or  in  military  applications  to 
assess  the  impact  of  the  earth's  atmosphere  on  military  systems.  For  these  kinds 
of  applications,  humidity  analyses  derived  in  the  manner  of  FNOC,  from  assimi¬ 
lation  systems  that  do  not  assimilate  humidity  observations,  can  develop  signifi- 

4  7 

cant  biases  owing  to  model  errors.  * 

Even  when  the  main  use  of  humidity  analysis  is  to  provide  initial  conditions 
for  a  large-scale  model,  the  assimilation  of  observed  humidity  data  is  important 
in  the  tropics.  There  the  coupling  of  the  mass  and  motion  fields  is  weak  and  the 
time  required  to  generate  reasonable  vertical  motion  fields  (and  hence  humidity 
fields)  is  several  days.  Furthermore,  improved  accuracy  in  describing  the  initial 
moisture  field  in  large-scale  models  has  been  shown  to  reduce  the  global  spin-up 
time  for  parameterized  processes  like  precipitation,  convective  heating,  and  sur- 
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face  moisture  fluxes.  For  these  reasons,  in  the  AFGL  global  modeling  effort, 
where  a  chief  goal  is  the  production  of  skillful  global  cloud  forecasts  in  the  one  to 
four  day  range,  humidity  analysis  and  assimilation  is  considered  a  necessary  and 
important  task. 
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In  this  study,  it  will  be  shown  that  the  choice  of  moisture  variable  in  the  ver¬ 
tical  interpolations  of  the  pre-  and  postprocessor  can  be  crucial  to  preserving  a 
reasonable  moisture  distribution  in  humidity  assimilations.  Section  2  reviews  the 
several  vertical  interpolation  methods  and  moisture  variable  choices  that  are  con¬ 
sidered.  Section  3  presents  results  of  the  use  of  relative  humidity  (RH),  dewpoint 
(T^),  and  specific  humidity  (Q)  in  applications  of  several  vertical  interpolations 
tested  within  the  NMC  pre-  and  postprocessor  applied  to  a  set  of  1978  FGGE  III-A 
analyses.  The  extent  to  which  the  vertical  interpolation  errors  of  the  pre-  and 
postprocessor  might  accumulate  in  a  4-D  assimilation  is  examined  in  Section  4  by 
means  of  a  "pseudo"  48-hr  data  assimilation.  Finally,  in  Section  5,  conclusions 
and  recommendations  are  presented. 


2.  MOISTURE  PREPROCESSING/POSTPROCESSING 

The  several  steps  in  a  given  forecast  cycle  in  Figure  1  are  depicted  in  further 
detail  in  Figure  2.  It  shows  that  the  preprocessing  step  entails  first  vertical  in¬ 
terpolation  (from  P  to  o)  and  then  spectral  transformation  from  physical  gridpoint 
space  to  spectral  coefficient  space  in  the  horizontal  direction.  The  postprocessor 
reverses  those  steps.  In  this  study,  primary  emphasis  is  given  to  the  errors  in 
the  vertical  interpolation  of  moisture  fields.  Some  attention  is  also  given  (in  Sec¬ 
tion  3.  3)  to  horizontal  truncation  errors  in  the  spectral  representation  of  moisture. 
This  study  begins  with  the  use  of  the  P  to  a  and  a  to  P  vertical  interpolation  algor¬ 
ithms  as  used  in  the  NMC  pre-  and  postprocessor.  This  choice  is  motivated  by 
the  inclusion  of  the  NMC  processors  in  the  family  of  new  NWP  models,  referred  to 
as  the  Advanced  Weather  Analysis  and  Prediction  System  or  AWAPS,  presently 
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being  implemented  at  AFGWC.  Concurrently,  several  research  requirements 
have  been  conveyed  to  AFGL  in  the  area  of  global  forecasting  and  analysis  of  mois¬ 
ture  and  clouds.  This  study  is  one  of  several  efforts  addressing  these  research 
requirements. 

The  general  configuration  of  the  global  4-D  assimilation  system  in  the  AWAPS 
is  given  in  the  top  of  Figure  1.  Through  an  interagency  agreement,  NMC  has  pro¬ 
vided  AFGWC  with  the  GSM,  NMI,  and  pre-  and  postprocessor.  In  addition 

AFGWC  has  developed  and  implemented  an  OI  objective  analysis  package  called  the 
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High  Resolution  Analysis  System  or  HIRAS,  Unlike  the  NMC  OI  analysis,  which 
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Figure  2.  Diagram  of  Components  of  a  Complete  Forecast  Cycle  (Solid-Line  Path) 
and  Components  of  a  Synthesized  Analysis  Sequence  (Dashed-Line  Path) 


analyzes  relative  humidity,  the  AFGWC  OI  analysis  model  is  presently  designed  to 
analyze  specific  humidity.  Consequently,  for  the  sake  of  expediency,  AFGWC  re¬ 
placed  relative  humidity  (RH)  with  specific  humidity  (Q)  in  the  vertical  interpola¬ 
tions  of  the  NMC  pre-  and  postprocessor. 

It  is  fairly  common  in  pre-  and  postprocessing  to  vertically  interpolate  the 
same  moisture  variable  that  is  objectively  analyzed.  This  is  apparent  in  Table  1, 
which  summarizes  for  five  NWP  facilities  the  moisture  variable  used  in  the  object¬ 
ive  analysis  and  vertical  interpolations  of  the  respective  4-D  assimilation  system. 
The  expedient  choice  of  moisture  variable  for  vertical  interpolation  is  not  neces¬ 
sarily  the  best  choice,  however.  Preliminary  AFGL  tests**  of  the  NMC/AFGWC 
pre  -  and  postprocessor  indicated  that  the  use  of  Q  in  the  vertical  interpolations  in¬ 
troduced  a  significant  positive  moisture  bias  in  the  6 -hr  first-guess  forecast  fields 
of  Figure  1.  This  problem  led  to  the  present  investigation  into  alternate  choices  of 
moisture  variable  and/or  vertical  interpolation  method  in  the  processing  steps. 


Table  1.  Summary  of  the  Analyzed  and  Vertically 
Interpolated  Moisture  Variables  in  the  4-D  Assim¬ 
ilation  at  Five  NWP  Facilities 


NWP 

Center 

Analyzed 

Moisture 

Variable 

V  ertically 
Interpolated 
Variable 

Vertical 

P-a/a-P 

Interpolation 

AFGWC8 

Q 

Q 

Linear 
in  LnP 

NMC2 

RH 

RH* 

Linear 
in  LnP 

ECMWF* 

Qt 

Q 

Cubic  spline 
in  LnP 

GFDL9 

Td 

Q 

Cubic  spline 
in  P 

GLAS10 

RH 

RH* 

Linear 
in  LnP 

’'‘Analyzed  residuals  of  RH  in  preprocessor 
^Layer  integrated 
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Except  for  ECMWF,  Table  1  is  current  at  the  time  of  writing.  In  the  case  of 

ECMWF,  Table  1  shows  the  procedures  used  to  perform  the  ECMWF  1979  FGGE 

ni-B  humidity  analyses.  As  suggested  in  Lorenc  and  Tibaldi,  *  the  latter  analyses 

likely  include  nontrivial  errors  from  vertical  interpolations  of  Q  (see  their  Figure 
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7).  As  documented  in  LOnnberg  and  Shaw,  ECMWF  has  since  replaced  cubic 
spline  interpolation  of  Q  in  LnP  with  a  procedure  using  layer  averages  of  RH  in  P. 
However,  the  present  authors  found  no  ECMWF  reports  in  the  public  domain  that 
documented  tests  leading  to  this  change.  The  present  study  will  show  why  vertical 
interpolation  of  Q  in  LnP  is  a  singularly  poor  choice  and  will  suggest  several  suit¬ 
able  alternatives. 

2.1  Vertical  Interpolation 

Vertical  interpolation  in  4-D  assimilation  is  necessary  whenever  the  analysis 

model  and  forecast  model  utilize  a  different  vertical  coordinate.  As  depicted  in 

Figure  1,  the  NMC  and  AFGWC  Ol  analyses  are  performed  on  constant  pressure 

surfaces  while  the  GSM  employs  the  standard  o-coordinate.  This  coordinate  is 

defined  by  a  *  P/P  ,  where  P  denotes  the  surface  pressure.  Hence  0  <o  s  1 
s  s 

and  P(a),  the  pressure  of  a  given  o-layer,  is  given  by  P(a)  =  aPs>  Table  2  pre¬ 
sents  the  o-layer  structure  used  in  the  NMC/AFGWC  GSM  for  various  vertical 
13 

resolutions.  Figure  3  schematically  illustrates  the  vertical  o-structure  of  the 
6-layer  GSM.  In  all  the  experiments  in  this  study,  we  employ  the  12-layer 
o-structure  of  Table  2,  following  the  NMC  operational  practice  at  the  time  re¬ 
ported  by  Sela.  ^ 

We  next  briefly  review  the  several  vertical  interpolations  to  be  tested  here  in 
the  NMC/AFGWC  pre-  and  postprocessor.  To  facilitate  concurrent  presentation 
of  these  vertical  interpolation  schemes,  a  generalized  notation  is  adopted.  One 
can  begin  by  considering  a  general  function  F  of  an  independent  variable  Z.  To 
first  order,  given  known  values  of  Fj  and  Fg  at  two  locations  and  Zg,  the  in¬ 
terpolated  value  F  at  some  location  Z  between  Z^  and  Zg  is  given  by  the  linear 
formula 

F  =  F1  +  A(Z  -  Zx)  ,  (1) 
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Figure  3.  Illustration  of  Vertical  o-Coordinate  in  NMC/AFGWC  Global 
Spectral  Model 


where 


A 


(2) 


The  error  or  remainder  term  implicit  in  Eq.  (1)  is 
(Z  -  Z  )2 

R  =  F"(Z*)  — - — —  .  (3) 

where  F"  is  the  second  derivative  of  F  at  some  Z*  between  Z  and  Z, . 

For  generality,  consider  F  and  Z  to  represent  elementary  transformations  of 
more  fundamental  physical  quantities  f  and  z,  such  that  F  =  F(f)  and  Z  =  Z(z).  In 
particular,  we  lef  f  generically  denote  the  chosen  moisture  variable  (RH,  Q,  or 
T^)  and  take  z  to  be  the  pressure  P.  In  Table  3,  it  is  seen  that  by  taking  F  and  Z 
to  represent  various  combinations  of  Lnf  or  f  and  LnP  or  P,  respectively,  the 
first-order  interpolation  algorithm  in  Eqs.  (1)  and  (2)  can  represent  any  one  of 
four  standard  interpolation  types  (namely  linear,  logarithmic,  exponential,  or 
power  law).  As  an  example,  the  equivalence  between  Eqs.  (1)  and  (2)  and  the  in¬ 
terpolation  type  given  in  Table  3,  in  the  case  of  power  law  interpolation,  is  shown 
in  Appendix  A. 


Table  3.  List  of  Interpolation  Type  Provided  by  Eqs.  (1)  and  (2)  for 
Various  Functional  Transformations,  F  =  F(f)  and  Z  =  Z(P),  of  the 
Moisture  Variable  f  and  Vertical  Coordinate  P 


Case 

F 

Z 

Interpolation  Type 

1 

f 

LnP 

logarithmic:  f  =  aLn(P/b) 

2 

f 

P 

linear:  f  =  a  P  +  b 

3 

Lnf 

LnP 

power  law:  f  =  a?*1 

4 

Lnf 

P 

exponential:  f  =aExp(P/b) 

In  the  NMC/AFGWC  pre-  and  postprocessor,  F  =  f  and  Z  =  LnP  (case  1  in 
Table  3).  Thus,  the  NMC/AFGWC  processors  assume  the  moisture  variable  f 
varies  logarithmically  with  P.  This  assumption  is  common  in  Table  1.  We  shall 
find  in  Section  3.  2  that  this  is  a  good  assumption  for  the  choice  of  f  =  Td  and  a  fair 
assumption  for  f  =  RH.  For  f  =  Q  on  the  other  hand,  we  shall  find  that  Q  does  not 
vary  logarithmically  with  P,  but  rather  closely  satisfies  a  power  law  in  P. 


The  special  case  of  extrapolation  must  also  be  considered.  In  this  case  levels 
(Zj,  Z2)  would  be  the  two  lowest  or  two  highest  levels  at  which  F  is  known,  and 
thus  either  Z<  Zj  <  Z^  (upper  extrapolation)  or  Zj  <  Z^  <  Z  (lower  extrapola¬ 
tion).  Different  extrapolation  assumptions  are  applied  to  the  moisture  variable  f 
in  the  NMC/AFGWC  preprocessor  and  postprocessor.  In  the  preprocessor,  linear 
extrapolation  following  Eqs.  (1)  and  (2)  is  applied,  while  in  the  postprocessor  con¬ 
stant  extrapolation  is  used.  Both  the  interpolation  and  extrapolation  procedures 
for  the  moisture  variable  f  in  the  NMC/AFGWC  processors  are  summarized  in 
Table  4  and  illustrated  schematically  in  Figure  4. 

Table  4.  Summary  of  Vertical  Interpolation  and  Extrapolation  Procedures  for  the 
Moisture  Field  in  the  NMC/AFGWC  Pre-  and  Postprocessor  (Where  Z  =  LnP) 


Preprocessor 

(P-a) 

Postprocessor 

(or-P) 

6  standard  pressure 
levels* 

GSM  moist  a  layers  t 

GSM  moist  a  layers  t 

6  standard  pressure 
levels* 

Linear  in  Z 

Eqs.  (1)  and  (2) 

Linear  in  Z 

Eqs.  (1)  and  (2) 

Linear  in  Z 

Eqs.  (1)  and  (2) 

Constant 

F(Z)  =  F(Zfa) 

Linear  in  Z 

Eqs.  (1)  and  (2) 

Constant 

F(Z)  =  F(Za) 

Levels  of  Known  Data 


Levels  of  Interpolated 
Data 

Interpolation 

Lower  Extrapolation 

Upper  Extrapolation 
2<  Za<  ^ 


noo.  85,  70.  50,  40,  30  cb 
tSee  Table  2. 

2.2  Moisture  Variables 

In  this  study,  relative  humidity  (RH),  specific  humidity  (Q),  and  dewpoint 
(T^)  are  tested  in  the  vertical  interpolations.  It  is  evident  in  Table  1  that  these 
moisture  variables  are  current  choices  in  global  objective  analysis  models  and 
thus  likely  candidates  for  vertical  interpolation. 

Fundamental  expressions  for  RH,  Q,  and  Td  are  derived  from  the  vapor 
pressure,  e,  and  saturation  vapor  pressure,  es<  given  with  respect  to  water  by 


es(T)  =  C/EXPI  L/(1.61RdT)] 
e(Td)  =  C/EXP[  L/(1.61RdTd)] 
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Figure  4.  Graphical  Illustration  of  the  Interpolation  and  Extrapolation  of 
Moisture  Variable  f  in  the  NMC/AFGWC  Pre-  and  Postprocessor  (Z  =  LnP) 
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Here  the  constant  C  =  2.  645  X  10  cb,  the  latent  heat  of  vaporization 
L  =  2.  51  X  10®  JKg  *  (at  T  =  273.  16K),  and  the  gas  constant  for  dry  air 
Rd  =  287  JKg  *K  *.  RH  and  Q  are  obtained  from  e  and  eg  through 

RH(Td,T)  =  (e/es)  X  100  .  (6) 

Q(Td'P)  =  P  -  O.^Se  '  (7) 

Lastly,  is  obtained  by  rewriting  Eq.  (5)  according  to 
L/(1.61R  ,) 

Td  2  Ln(C/e)  *  (8) 

With  Eqs.  (4)  through  (8)  and  given  a  value  of  T  and  P,  any  value  of  RH,  Q,  or 
can  be  used  to  derive  either  of  the  remaining  two  moisture  variables  whenever 
needed.  While  polynomial  approximations*8  for  e  and  e  can  be  used  in  place  of 
Eqs.  (4)  and  (5),  we  chose  not  to  pursue  these  approximations  here.  Implicit  in 
Eqs.  (4)  through  (8)  is  the  computation  of  the  humidity  variables  with  respect  to 
vapor  pressure  over  water  for  all  values  of  temperature,  regardless  of  whether 
or  not  it  is  subfreezing.  This  is  necessary  to  be  consistent  with  the  FGGE  III -A 
objective  analyses  of  RH,  which  are  specified  universally  with  respect  to  water. 

As  discussed  below,  they  are  used  here  for  both  input  and  assessment  fields. 

2.3  The  Processing  Sequence 

In  a  complete  4-D  assimilation  system  such  as  that  in  Figure  1,  the  errors 
arising  from  the  pre-  and  postprocessing  steps  can  be  difficult  to  separate  from 
other  sources  of  error.  Therefore,  to  highlight  moisture  processing  error.  Sec¬ 
tion  3  examines  the  processing  error  obtained  by  generating  and  verifying  a  "syn¬ 
thesized  analysis".  As  depicted  by  the  dashed-line  path  in  Figure  2,  the  synthe¬ 
sized  analysis  is  produced  by  applying  the  NMC/AFGWC  pre-  and  postprocessor 
alone  to  an  input  FGGE  III-A  analysis.  As  such,  the  synthesized  analysis  repre¬ 
sents  a  static  transformation  of  the  original  input  analysis,  since  no  forecast 
model  (or  initialization)  is  executed. 

Table  5  gives  the  specific  steps  in  the  cycle  of  pre-  and  postprocessing  used 
in  this  study  to  generate  a  synthesized  analysis.  We  note  that  30-wave  rhomboidal 
truncation  (denoted  R30)  is  utilized  in  the  spectral  transform  steps  of  Table  5, 


15.  Lowe,  P.  (1977)  Approximating  polynomial  for  the  computation  of  saturation 
vapor  pressure,  J.  Appl.  Meteorol.  16:100-105. 
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Table  5.  Sequence  of  Moisture  Computations  in  the  NMC/AFGWC 
Preprocessor  (a)  and  (b)  Postprocessor 


(a)  Preprocessor 

Steo 

Computation  Performed 

1 

Input  FGGE  III-A  analysis  of  RH  and  heights  H  on  2. 5°  X  2.5°  (144  X  73) 
latitude -longitude  grid  on  standard  pressure  surfaces. 

2 

Horizontally  interpolate  the  analyzed  RH  and  H  fields  to  the  76  GSM 
Gaussian  latitudes  on  standard  pressure  surfaces. 

3 

Hydrostatically  derive  temperature  T  (and  surface  pressure  Ps)  on 
standard  pressure  surfaces  from  the  height  fields,  H,  and  input  model 
terrain. 

4 

Apply  Eqs.  (4)  through  (8)  as  necessary  to  convert  the  (144  X  76)  fields 
of  T  and  RH  to  the  desired  moisture  field  f  =  (RH,  Q,  or  Tj)  on  standard 
pressure  surfaces. 

5 

Vertically  interpolate  moisture  field  f  of  step  4  from  standard  pres¬ 
sure  surfaces  to  GSM  ct -layers  of  Table  2  using  one  of  the  linear 
interpolation/extrapolation  types  in  Table  3.  Hydrostatically  derive 

T  on  a-layers  using  T  and  H  on  standard  pressure  surfaces. 

6 

Apply  Eqs.  (4)  through  (8)  as  necessary  to  derive  Q  from  T  and  moisture 
field  f  =  (RH,  Q,  or  T^)  on  GSM  o-layers. 

7 

Spectrally  transform  (144  X  76)  gridded  fields  of  Q  and  T  on  GSM 
a-layers  to  spherical  harmonic  coefficients  of  Q  and  T. 

(b)  Postprocessor 

Step 

Computation  Performed 

1 

Spectrally  transform  spherical  harmonic  coefficients  of  Q  and  T  to 
(144  X  76)  Gaussian  latitude  grid  on  GSM  a-layers. 

2 

Apply  Eqs.  (4)  through  (8)  as  necessary  to  convert  gridded  Q  and  T 
fields  on  GSM  a-layers  to  moisture  field  f  =  (RH,  Q,  or  Tj).  Choice 
of  f  here  agrees  with  that  in  step  4,  Table  5(a). 

3 

Vertically  interpolate  T  and  moisture  field  f  from  GSM  a-layers  to 
standard  pressure  surfaces  using  the  linear  interpolation/extrapolation 
type  of  Table  3  used  in  step  5,  Table  5(a). 

4 

Apply  Eqs.  (4)  through  (8)  as  necessary  to  convert  T  and  moisture  field 
f  to  (144  X  76)  RH  field  on  standard  pressure  surfaces.  (This  step  un¬ 
necessary  if  moisture  field  f  =  RH. ) 

5 

Output  (144  X  76)  RH  field  on  Gaussian  latitude  grid  on  standard  pressure 
surfaces  for  subsequent  verification. 
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again  following  the  NMC  operational  practice  at  the  time  reported  by  Sela.  In 
some  tests  in  Section  3.3,  the  spectral  transforms  in  Table  5  are  eliminated  to 
isolate  processing  errors  due  to  vertical  interpolation. 


3.  ANALYSIS  OF  PROCESSING  ERROR 


In  Figure  2,  the  solid-line  path  depicts  the  total  forecast  cycle  F^.  The  total 
forecast  error  shall  be  denoted  (F^,  A^)  representing  the  verification  of  the  total 
forecast  cycle  F1  against  the  analysis  A^,  both  valid  at  time  tj.  The  dashed-line 
path  in  Figure  2  depicts  the  synthesis  or  static  processing  cycle  S^.  The  synthe¬ 
sis  or  processing  error  shall  be  denote'*  (S^,  A^)  representing  the  verification  of 
the  synthesized  analysis  against  the  original  analysis  A^.  This  section  will  ex¬ 
amine  the  processing  error  (S^,  A^),  while  Section  4  will  consider  total  forecast 
error  as  well. 

In  this  study  the  two  error  quantities  computed  at  each  standard  pressure 
level  are  the  global  mean  RMSE  and  BIAS.  They  are  obtained  from 
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where  are  the  "model-produced"  RH  values  (either  F^  or  S^)  and  are  the 
"verifying"  RH  values  of  the  analysis  A^  at  each  horizontal  grid  point  (i,  j)  on  a 
given  pressure  surface  P.  The  indices  i=l,  2,  ....  I  and  J=l,  2,  ...»  J  respect¬ 
ively  denote  the  longitudinal  points  (equally  spaced)  and  the  latitudinal  points  (non- 

equally  spaced  at  Gaussian  latitudes)  of  the  global  grid.  W.  is  the  Gaussian  quad- 

J  3 

rature  weight  satisfying  £  =  2.  RMSE.  and  BIAS,  represent  zonal  mean 

j=l  J  13 

values. 

The  1978  FGGE  III-A  analyses,  which  are  products  of  the  operational  NMC 
4-D  data  assimilation  system  in  1978,  were  used  as  the  input  and  assessment  ob¬ 
jective  analyses  for  this  study.  These  analyses  provide  fields  of  height  H  at  the 
lowest  12  standard  pressure  levels  and  RH  at  the  lowest  six  standard  pressure 


levels  on  a  regular  144  X  73  latitude -longitude  grid  (2.  5°  X  2.  5°).  The  III-A 
analyses  were  interpolated  horizontally  to  the  76  Gaussian  latitudes  of  the  R30 
GSM,  transforming  the  fields  to  144  X  76.  Throughout  our  verifications,  we  use 
the  RH  analyses  [step  2,  Table  5(a)]  and  RH  forecasts  [step  5,  Table  5(bl]  on 
this  144  X  76  Gaussian  grid.  Thus  in  Eqs.  (9)  and  (10),  I  =  144  and  J  =  76.  By 
staying  consistently  on  the  Gaussian  grid,  we  eliminate  horizontal  interpolation 
error  from  our  measures  of  vertical  interpolation  error. 

Note  that  the  quantities  M.^  and  V..  in  Eqs.  (9)  and  (10)  shall  always  be  taken 
as  values  of  RH,  regardless  of  the  moisture  variable  used  in  the  vertical  interpo¬ 
lations.  Thus,  the  vertical  interpolations  of  Q  and  T^,  when  verified  in  terms  of 
RH,  may  suffer  additionally  from  processing  errors  in  the  temperature  field 
[  step  4,  Table  5(b)].  While  doing  so  may  risk  representing  the  performance  of 
RH  interpolations  too  favorably,  it  is  easier  to  judge  the  significance  of  moisture 
field  errors  when  expressed  in  terms  of  RH.  As  will  be  seen  later,  the  implicit 
favoritism  toward  RH-based  interpolations  in  the  verifications  will  prove  to  be 
small  in  most  cases  and  never  a  dominant  factor  in  the  comparisons.  Neverthe¬ 
less,  for  completeness,  at  the  end  of  Section  3  we  try  to  quantify  the  degree  of 
RH -favoritism  in  the  verifications. 

3. 1  Extrapolation  Error 

This  section  begins  the  examination  of  verification  statistics  for  examples  of 
processing  error  (S^  A^.  The  first  examples  show  substantial  extrapolation 
errors  at  the  upper  and  lower  levels,  so  we  initially  focus  on  these  errors  and 
methods  to  eliminate  them.  The  investigation  of  extrapolation  errors  in  this  sec¬ 
tion  will  entail  only  the  standard  logarithmic  interpolation  type  (case  1,  Table  3). 
The  next  section,  which  examines  interpolation  errors  in  the  interior  of  the  moist 
vertical  domain,  shall  pursue  the  remaining  interpolation  types  in  Table  3. 

It  is  important  for  this  section  to  note  in  Table  2  that  NMC  presently  does  not 
carry  a  prognostic  moisture  field  at  all  a-layers.  Rather,  to  save  computation 
time,  NMC  chooses  in  its  GSM  to  carry  a  prognostic  moisture  field  only  at  cr-layers 
essentially  in  the  troposphere.  It  is  evident  in  Table  2  that  for  typical  surface 
pressures,  the  highest  moist  a-layer  that  NMC  recommends  is  neither  consistently 
above  nor  below  30  cb  for  all  three  vertical  resolutions.  For  example,  in  the  9- 
and  12 -layer  resolutions  (in  which  the  first  six  cr-layers  are  identical)  six  and 
seven  moist  a-layers  are  recommended  respectively.  Therefore,  to  further  save 
computation  time  in  the  GSM,  we  first  tested  using  just  six  moist  a-layers  in  our 
12 -layer  experiments  with  the  pre-  and  postprocessor. 

For  the  three  choices  of  RH,  Q,  or  T^  as  the  vertically  interpolated  moisture 
variable  f.  Figure  5  presents  the  vertical  profiles  of  the  global  mean  BIAS  and 
RMSE  for  the  processing  error  (Sj,  A^)  obtained  by  applying  the  processing  steps 
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Figure  5.  Vertical  Profiles  of  Global  Mean  BIAS  and  RMSE  for  the  (Si,  A^) 
Verification  of  the  00Z  FGGE  III -A  Analysis  of  17  January  1978,  Synthesized 
by  Applying  the  Computations  of  Tables  4  and  5  to  Six  Moist  a-Layers  for  the 
Three  Choices  of  Vertically  Interpolated  Moisture  Variable  of  F  =  (RH,  Q,  T^) 


of  Table  5  to  the  00Z  FGGE  III -A  analysis  of  17  January  197  8.  The  results  in 
Figure  5  incorporate  the  operational  NMC/AFGWC  first-order  interpolation/ ex¬ 
trapolation  algorithms  of  Eqs.  (1)  and  (2)  using  F  =  f  and  Z  =  LnP  (case  1,  Table  3), 
as  further  depicted  in  Table  4  and  Figure  4. 

In  Figure  5,  the  strikingly  large  BIAS  errors  at  100  and  30  cb  in  the  Q-  and 
T^ -based  profiles  result  from  a  common  extrapolation  problem,  which  we  illus¬ 
trate  shortly.  (The  reader  should  note  the  break  in  horizontal  scale  at  the  right 
of  each  plot  in  Figure  5. )  Though  much  smaller,  the  bias  error  at  100  and  30  cb 
in  the  RH-based  profile  stems  from  the  same  extrapolation  weakness.  Inspection 
of  the  latitude -height  cross-sections  (not  shown)  of  the  zonal  mean  BIAS  corres¬ 
ponding  to  Figure  5  confirmed  that  both  the  100-  and  30-cb  BIAS  were  quite  uniform 
across  all  latitudes  and  not  dominated  by  extreme  values  in  some  limited  latitude 
zone. 

The  left  side  of  Figure  6  schematically  illustrates  that,  given  the  input  of 
moisture  values  on  only  six  standard  pressure  surfaces  and  given  a  moisture  var¬ 
iable  that  decreases  with  height  in  the  upper  troposphere,  then  the  constant  upper 
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Figure  6.  Illustration  of  Bias  Effects  of  Currently  Operational 
Upper  Moisture  Extrapolation  Sequence  in  NMC/AFGWC  Pre- 
and  Postprocessor  for  6  or  7  Moist  o-Layers  in  Either  the  9- 
or  12 -Layer  Resolutions  in  Table  2.  The  moisture  variable  f 
is  assumed  to  decrease  with  height  (see  Figure  7) 

extrapolation  in  the  postprocessor  will  always  yield  a  positive  (S1,  A1)  bias  at 
30  cb  if  the  highest  moist  o-layer  is  below  30  cb.  This  conclusion  holds  regard¬ 
less  of  the  total  number  of  moist  o-layers  in  the  chosen  vertical  resolution.  The 
magnitude  of  this  30-cb  bias  will  depend  on  the  vertical  "lapse  rate"  of  the  chosen 
moisture  variable,  the  proximity  of  the  highest  moist  o-layer  to  30  cb,  and  whether 
the  moisture  variable  is  a  relative  or  absolute  measure  of  humidity.  Figure  7 
presents  the  sample  mean  vertical  profile  of  the  globally  averaged  RH,  Q,  and 
for  six  January  1978  FGGE  UI-A  analyses  (original,  nonsynthesized).  For  refer¬ 
ence  purposes.  Figure  7  also  gives  the  sample  mean  temperature  profile,  as  well 
as  representative  locations  of  the  six  moist  a-layers  (for  Pg  =  101.3  cb).  We  find 
that  all  three  moisture  variables  decrease  with  height  at  all  levels.  Similarly, 
based  on  the  mean  temperature  profile,  the  mean  profiles  of  saturation  Q  and  T^ 
(not  shown)  decrease  with  height  at  all  levels.  Obviously,  saturation  values  of  Q 
become  very  small  at  the  cold  temperatures  at  30  cb,  rendering  interpolation  or 
extrapolation  errors  in  Q  (when  expressed  as  RH  errors)  very  sensitive  to  small 
errors.  For  the  current  o-structure  (12-layer)  and  the  "climatology"  of  Figure  7, 
the  positive  (S^,  A^)  bias  at  30  cb  in  Figure  5  is  acceptably  small  when  using  RH 
(NMC  choice),  but  is  unacceptably  large  when  using  Q  (AFGWC  choice)  or  T^. 

These  large  biases  are  shown  in  Figure  5  as  equivalent  RH  biases. 
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Figure  7.  Vertical  Profiles  of  Global  Mean  T,  Tj,  RH,  and  Q  for  a 
Sample  of  Six  Winter  1978  FGGE  III-A  Analyses  for  00Z  on  15,  17, 

19  January  and  12 Z  on  15,  18,  20  January.  Location  of  lowest  six 
a-layers  of  9-  or  12-layer  resolution  in  Table  2  is  shown  for 

P  =  101,3  cb 
s 

To  eliminate  the  upper  extrapolation  bias,  the  approach  illustrated  on  the 
right  of  Figure  6  is  recommended,  wherein  the  highest  moist  o-layer  is  chosen 
above  30  cb.  Thereby,  upper  extrapolation  of  moisture  is  invoked  only  in  pre¬ 
processing  (which  employs  linear  instead  of  constant  extrapolation.  Figure  4) 
while  in  the  postprocessor  the  final  30 -cb  moisture  value  is  now  determined  by 
interpolation.  In  fact,  given  a  moist  o-layer  above  30  cb  and  at  least  one  between 
30  through  40  cb,  the  processing  sequence  on  the  right  of  Figure  6  will  exactly 


reproduce  the  original  input  moisture  value  f  at  30  cb  (disregarding  roundoff  error 
and  horizontal  spectral  truncation  error). 

Rather  than  follow  the  right  of  Figure  6,  one  may  argue  that  instead  one  could 
still  retain  only  six  moist  a  -layers  and  employ  linear  instead  of  constant  upper 
extrapolation  in  the  postprocessor.  This  is  a  fair  alternative  for  RH  or  Td  (though 
not  Q),  as  shown  shortly.  But  this  alternative  does  not  exactly  reproduce  the  input 
moisture  value  at  30  cb  unless  at  least  two  moist  cr-layers  exist  between  30  through 
40  cb.  Finally,  we  caution  that  in  coarse  vertical  resolutions,  where  the  first 
a-layer  above  30  cb  may  exist  well  above  the  tropopause  (above  which  the  moisture 
content  in  the  stratosphere  is  thought  to  become  quasi-constant  with  height)  or  in 
resolutions  carrying  several  active  moist  a-layers  in  the  stratosphere,  other  ex¬ 
trapolation  approaches  to  replace  that  recommended  in  the  right  of  Figure  6  should 
be  explored. 

The  100-cb  negative  bias  in  all  three  profiles  in  Figure  5  arises  in  similar 
fashion  from  the  constant  lower  extrapolation  in  postprocessing  shown  in  Table  4 
and  Figure  4.  The  left  of  Figure  8  illustrates  the  drying  effect  of  constant  lower 
extrapolation  of  a  moisture  variable  that  decreases  with  height.  This  drying  effect 
is  especially  pronounced  in  Figure  5  for  the  absolute  moisture  variables  of  f  =  Q 
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Figure  8.  Illustration  of  Bias  Effects  of  Current  (Left)  and  Pro¬ 
posed  (Right)  Lower  Moisture  Extrapolation  Sequence  in  NMC/ 
AFGWC  Pre-  and  Postprocessor  in  Either  the  9-  or  12 -Layer 
Resolutions  in  Table  2.  The  moisture  variable  f  is  assumed  to 
decrease  with  height  (Figure  7) 


or  £  =  T^,  since  the  postprocessor  applies  linear  lower  extrapolation  of  tempera¬ 
ture  and  thus  the  saturation  values  of  Td  and  Q  will  increase  downward.  The  right 
of  Figure  8  shows  the  improvement  obtained  by  applying  linear  lower  extrapolation 
for  the  moisture  variable  in  postprocessing.  Indeed,  whenever  two  or  more  o- 
layers  are  carried  within  85  to  100  cb,  linear  lower  extrapolation  in  postprocess¬ 
ing  will  give  100-cb  synthesized  moisture  values  f  that  agree  exactly  with  the  input 
values  (again  disregarding  roundoff  and  spectral  truncation  errors). 

Figure  9  shows  the  results  of  repeating  the  test  in  Figure  5  using  the  recom¬ 
mended  linear  lower  and  upper  extrapolation  of  the  moisture  variable  throughout 
the  preprocessor  and  postprocessor.  The  solid  curves  of  Figure  9  were  obtained 
using  seven  moist  a-layers  (following  the  right  of  Figure  6),  while  the  dashed 
curves  were  obtained  using  six  moist  a-layers.  In  the  case  of  seven  moist  layers, 
the  BIAS  at  100  and  30  cb  is  now  virtually  zero  for  all  the  moisture  variables.  In 


the  case  of  six  moist  layers,  the  BIAS  and  RMSE  results  are  acceptably  compar¬ 
able  to  those  of  the  case  of  seven  moist  layers,  except  in  the  Q-based  vertical  in¬ 
terpolations.  In  the  latter  case  the  nonlinearity  of  the  sample  mean  Q-profile  with 
respect  to  LnP  in  the  upper  troposphere  in  Figure  7  shows  Q  to  be  less  amenable 
than  RH  or  to  linear  upper  extrapolation  in  LnP.  Figure  10  repeats  the  test  of 
Figure  9,  with  similar  results,  for  the  00Z  FGGE  III-A  analysis  of  18  July  1978. 
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Figure  10.  As  in  Figure  9  Except  for  the  00Z  FGGE  III-A  Analysis  of  18  July  1978 


From  the  results  of  this  section,  it  is  clear  that  the  upper  and  lower  moisture 
extrapolation  in  the  current  NMC/AFGWC  processors  is  acceptable  (though  not 
optimal)  for  f  =  RH,  but  is  clearly  unacceptable  for  f  *  Q  or  f  =  (Figure  5). 
Admittedly,  in  the  original  design  of  the  NMC  processors  the  intent  may  have  been 
to  apply  the  moisture  interpolation/extrapolation  algorithms  only  to  f  =  RH.  Un¬ 
fortunately,  this  likely  intent  was  not  followed  in  the  later  adaptation  of  the  NMC 
processors  for  use  by  AFGWC. 

In  conclusion,  for  any  choice  of  moisture  variable  f  =  (RH,  Q,  or  T^)  in  the 
NMC/AFGWC  processors,  it  is  recommended  that  (1)  linear  rather  than  constant 
upper  and  lower  extrapolation  be  applied  in  postprocessing  and  (2)  the  highest 


moist  o-layer  be  the  lowest  layer  consistently  above  30  cb.  In  Table  2  then  the  use 
of  5,  7,  and  7  moist  o-layers  is  recommended  in  the  6-,  9-,  and  12-layer  resolu¬ 
tions,  respectively.  The  second  recommendation  above  is  not  crucial  for  f  =  RH 
or  T^,  but  is  imperative  for  logarithmic  interpolations  using  f  =  Q.  Lastly, 
results  later  in  the  next  section  suggest  that  six  moist  o-layers  is  sufficient  also 
when  using  f  =  Q  provided  F  =  Lnf  =  LnQ  and  Z  =  LnP  (power  law  interpolation), 
since  LnQ  will  be  shown  to  be  nearly  linear  in  LnP. 

3.2  Interpolation  Error 

In  the  previous  subsection,  significant  reductions  in  vertical  extrapolation 
error  were  achieved.  Yet,  a  further  reduction  in  the  processing  error  still  pres¬ 
ent  in  the  solid  curves  of  Figures  9  and  10  is  desirable.  In  particular,  one  would 
like  to  reduce  the  distressingly  large  positive  BIAS  and  RMSE  in  the  40  to  70  cb 
domain  of  the  Q-based  error  profiles.  Toward  this  end,  this  section  compares 
the  performance  of  the  standard  logarithmic  interpolation  (F  =  f  and  Z  =  LnP)  with 
that  of  the  other  interpolation  types  in  Table  3. 

We  recall  from  Eq.  (3)  that  to  achieve  accuracy  from  the  first-order  inter¬ 
polation  of  Eqs.  (1)  and  (2),  the  dependent  variable  F  must  be  fairly  linear  in  Z 
(that  is,  F"  *  0).  If  F(Z)  is  not  linear  but  rather  consistently  convex  (concave), 
as  given  by  F"  >  0  (F"  <  0)  throughout  the  domain,  then  a  successive  pair  of  first- 
order  interpolations  from  P  -*  a  and  a  -*  P  will  yield  a  positive  (negative)  bias. 

This  is  demonstrated  schematically  in  Figure  11.  It  is  clear  from  Figure  7  that 
for  F  =  f  =  Q  and  Z  =  LnP  (logarithmic  interpolation),  the  January  mean  function 
F(Z)  is  consistently  convex.  It  will  be  shown  that  the  convex  nature  of  Q  in  Z  =  LnP 
is  the  direct  cause  of  the  positive  bias  at  the  interior  pressure  levels  in  the 
Q-based  profile  of  Figures  9  and  10. 

One  can  anticipate  the  preferable  interpolation  type  for  Q-based  interpolations 

by  examining  Figure  12,  which  shows  the  N.  H.  1978  FGGE  III-A  January  and  July 

sample  mean  vertical  profiles  of  Q  for  all  four  specifications  of  F  and  Z  in  Table  3. 

Figure  13  shows  the  corresponding  observed  mean  Q-profiles  obtained  from  the 

1  fi 

10-year  (1963-1973)  climatological  RAOB  statistics  of  Oort.  In  an  absolute 
sense,  the  FGGE  III-A  profiles  of  Q  are  too  moist  compared  to  the  RAOB  clima¬ 
tology  (except  at  100  cb  in  January),  but  the  shapes  of  the  vertical  profiles  in 
Figures  12  and  13  are  in  good  agreement.  (In  the  case  of  LnQ,  Figures  12  and  13 
show  LnQ;  that  is,  Ln  of  the  hemispheric  average  of  Q.  We  would  prefer  to  show 
LnQ,  the  hemispheric  average  of  LnQ,  but  this  was  not  available  in  the  Oort 


16.  Oort,  A.  (1983)  Global  Atmospheric  Circulation  Statistics,  1958-1973,  NOAA 

S.  Government  Printing  Office.  Washington. 
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Figure  11,  Illustration  of  Relation  of  Convex  and  Concave  Curvature 
of  a  General  Monotonic  Increasing  Function  F(Z)  to  Bias  in  First- 
Order  Vertical  Interpolation  in  Pre-  and  Postprocessing 


statistics.  Strictly  speaking  LnQ  *  LnQ,  but  direct  calculations  of  LnQ  profiles 
from  the  FGGE  III -A  data  showed  a  nearly  uniform  displacement  of  LnQ  from  the 
profile  of  LnQ  at  all  levels.  Therefore,  LnQ  preserves  the  essential  shape  or 
curvature  of  LnQ. ) 

For  either  season.  Figures  12  and  13  show  that  on  the  average  Q  is  sharply 

convex  in  Z  =  LnP  and  modestly  convex  in  Z  =  P,  while  LnQ  is  nearly  linear  in 

Z  =  LnP  and  modestly  concave  in  Z  =  P.  Thus,  in  vertical  interpolations  of  Q 

using  Eqs.  (1)  and  (2),  we  expect  from  Figures  11  through  13  that  the  logarithmic, 

linear,  power  law,  and  exponential  interpolations  of  Table  3,  respectively,  will 

yield  biases  that  are  positive,  smaller  positive,  relatively  zero,  and  negative.  In 

conjunction  with  Figure  13,  a  least  squares  fit  of  the  six  discrete  (Q, P)  values  in 

Figure  13(a)  was  performed  against  a  logarithmic,  linear,  exponential,  and  power 

law  function  in  turn.  These  curve  fitting  tests  yielded,  respectively,  squared 

correlation  coefficients  of  0.  735,  0.851,  0.978,  and  0.  998.  Hence,  the  intrinsic 

mean  vertical  variation  of  Q  in  P  is  best  described  by  a  power  law,  which  in  the 

3  5 

above  exercise  was  given  by  Q  =  aP  '  . 

To  quantify  the  expected  response  in  bias  error  to  variations  in  the  curvature 
of  F(Z)  for  Q-based  interpolations,  an  idealized  test  was  performed  that  provided 
Table  6(a).  In  Table  6(a)  are  the  results  of  applying  Eqs.  (1)  and  (2),  including 
first -order  extrapolation  whenever  necessary,  for  each  case  in  Table  3  to  a  pair 
of  P  —  o  and  o  —  P  interpolations  of  the  single  January  hemispheric  mean  profile 
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Figure  12.  _Vertical  Profiles  of  the  Analyzed  N.  H.  Mean  Specific 
Humidity,  Q.  and  LnQ,  Plotted  Against  P  and  LnP  for  an  Aver¬ 
aged  Sample  of  Six  1978  FGGE  III-A  Analyses  in  (a)  January 
(Using  00Z  on  15,  17,  19  January  and  12Z  on  15,  18,  20  January) 
and  (b)  July  (Using  00Z  on  16,  18,  20  July  and  12Z  on  17,  19,  22 
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of  Q  in  Figure  13(a).  In  these  interpolations,  the  seven  moist  cr-layers  of  the  12- 

layer  NMC  GSM  were  still  used,  with  P  =  101.  3  cb.  Using  the  January  hemi- 

s 

spheric  mean  temperature  profile  (noninterpolated)  corresponding  to  Figure  13(a), 
plus  Eqs.  (4)  through  (8)  as  necessary,  one  can  express  the  interpolation  error 
[Q(P)  -  Q(P)]  equivalently  as  an  RH  error,  as  given  in  Table  6(a). 

The  RH  errors  in  Table  6(a)  show  that  for  the  given  input  Q -profile,  the  sign 
and  relative  magnitude  of  the  interpolation  error  for  each  interpolation  type  is  as 
expected  from  Figure  11.  Most  notably,  the  standard  logarithmic  interpolation 
indeed  yields  the  largest  error,  while  power  law  interpolation  yields  the  smallest 
error.  Furthermore,  the  sign  of  the  error  changes  from  positive  to  negative  from 
left  to  right  in  Table  6(a)  in  agreement  with  the  change  from  convex  to  concave 
curvature  from  left  to  right  in  Figure  13(a). 

Even  the  cubic  spline  interpolation  scheme  does  not  eliminate  the  large  posi¬ 
tive  bias  arising  from  the  linear  interpolation  of  Q  in  LnP.  This  is  shown  in 

Table  6(b),  which  repeats  the  test  of  Table  6(a)  but  using  the  cubic  spline  interpo- 

17 

lation  scheme  described  in  Ahlberg  et  al  in  place  of  the  first-order  scheme  of 
Eqs.  (1)  and  (2).  Except  for  the  large  extrapolation  errors  at  30  cb,  which  we 
ignore  for  the  moment.  Table  6(b)  shows  that  for  each  functional  form  F(Z),  the 
cubic  spline  scheme  yields  RH  errors  only  negligibly  smaller  than  the  errors  of 
the  first-order  scheme  in  Table  6(a).  Of  particular  note,  the  first-order  inter¬ 
polation  of  LnQ  in  LnP  is  far  superior  to  cubic  spline  interpolation  of  Q  in  LnP. 
Thus,  Table  6  shows  that  by  choosing  an  appropriate  transformation  (f,  P)—(F,  Z) 
that  minimizes  F"  in  Eq.  (3),  first-order  interpolation  can  be  substantially  better 
than  cubic  spline  interpolation  applied  to  the  original  variables  (f, P).  Tables  6(c) 
and  6(d)  show  similar  results  and  correspond  to  Tables  6(a)  and  6(b),  respectively, 
but  for  the  July  mean  Q-profiles  in  Figure  13(b).  Finally,  with  the  same  implica¬ 
tions,  Table  7  repeats  the  tests  of  Table  6,  but  for  the  FGGE  III-A  January  and 
July  mean  Q-profiles  in  Figure  12.  Based  on  Tables  6  and  7,  the  use  of  cubic 
splines  was  not  pursued  further.  If  we  had  pursued  cubic  splines  we  presumably 
could  have  identified  methods  to  reduce  the  large  30  cb  errors  of  the  cubic  spline 

extrapolation,  which  in  Tables  6  and  7  utilized  "Type -3"  of  the  three  types  of  ex- 

17 

trapolation  or  end -conditions  that  Ahlberg  et  al  present  as  choices  with  the  cubic 
spline. 

It  is  instructive  to  consider  further  the  nature  of  first-order  interpolation 
errors.  It  holds  that  for  a  given  input  profile  F(Z)  of  a  given  curvature  (concave 
or  convex),  the  application  of  n  repetitions  of  a  pair  of  P  —a  and  ct  — P  interpola¬ 
tions  using  Eqs.  (1)  and  (2)  will  yield  an  interpolated  profile  that  converges  to  a 


17.  Ahlberg,  J.  ,  Nilson,  E. ,  and  Walsh,  J.  (19S7)  The  Theory  of  Splines  and 
Their  Applications,  Academic  Press,  New  York,  284  pp. 
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Table  7.  As  in  Table  6  Except  FGGE  1II-A  Analyzed  N.H.  Mean  Q-Profiles  of  Figure  12  Replace  Those  of  Figure 
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linear  profile  (F  -*  0)  for  large  n.  Figure  11  showed  the  mechanics  of  this  effect 
and  Figure  14  provides  an  example.  In  Figure  14,  the  F(Z)  profile  of  F  =  Q, 

Z  =  LnP  in  Figure  13(a)  was  subject  to  n  =  8,  16,  and  24  iterations  of  the  P  —  a  and 
o-»  P  logarithmic  interpolation  (corresponding  respectively  to  the  number  of  pre- 
and  postprocessor  pairs  executed  in  a  2-,  4-,  and  6-day  run  of  the  assimilation 
system  of  Figure  1).  In  actual  practice,  the  asymptotically  linear  profile  in  Fig¬ 
ure  14  may  not  be  reached  owing  to  checks  for  supersaturation.  Nevertheless,  to 
inhibit  climatic  drift  in  humidity  assimilations  due  to  vertical  interpolations.  Fig¬ 
ure  6  clearly  demonstrates  that  one  should  choose  a  moisture  function  F(Z)  that  on 


Q(g/kg)  LINEAR  SCALE 

Figure  14.  A  Sample  of  Interpolated  Q-Profiles  for  n  Iterations 
(n  =  0,  8,  16,  24)  of  the  Pair  of  P  -*  o  and  o-  P  Logarithmic  In¬ 
terpolations  [Eqs.  (1)  and  (2)  With  F  =  Q,  Z  =  LnP)  Applied  to 
the  Input  Q-Profile  of  Table  6(a)  and  Using  the  Seven  Moist  a- 
Layers  of  the  12 -Layer  Column  of  Table  2  With  Ps  =  101.  3  cb 


The  interpolation  tests  in  Tables  6  and  7  were  idealizations,  since  each  was 
performed  on  a  single  hemispheric  mean  profile  of  Q.  Certainly,  local  vertical 
profiles  of  Q  within  a  global  moisture  analysis  are  less  smooth,  possessing  more 
vertical  variability  and  changes  in  curvature  than  their  hemispheric  mean  counter' 
parts.  Therefore,  to  establish  the  representativeness  of  Tables  6  and  7  for  prac¬ 
tical  applications,  we  extended  the  tests  of  the  four  interpolation  types  to  actual 
executions  of  the  pre-  and  postprocessor  on  a  global  data  set  and  to  all  three 
moisture  variables.  Table  8  shows  the  (Sj,  A^  BIAS  and  RMSE  obtained  for  each 
choice  of  f  =  (RH,  Q,  or  Td)  and  interpolation  type  in  Table  3  used  within  the  pre- 
and  postprocessor  executed  on  the  same  input  FGGE  analysis  as  in  Figure  9.  For 


Table  8.  Vertical  Profiles  of  the  Global  Mean  BIAS  (a)  and  RMSE  (b)  for  the  Synthesized  FGGE  III-A  Analysis 
(17  January  1978,  00Z)  Obtained  Using  the  Four  Interpolation  Types  of  Table  3  and  the  Moisture  Variables 
f  =  (RH,  Q,  T<j)  in  the  Vertical  Interpolations  of  the  NMC/AFGWC  Pre-  and  Postprocessor 
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reference,  the  leftmost  error  column  in  Table  8(a)  and  (b)  corresponds  to  the 
solid  curves  in  Figure  9. 

Overall,  Table  8  shows  that  the  Q-based  processing  error  is  more  sensitive 
to  interpolation  type  than  that  of  RH  or  T^.  Notably,  the  magnitude  and  sign  of 
the  Q-based  errors  follow  closely  those  obtained  earlier  in  the  idealized  tests  of 
Tables  6  and  7  --  showing  a  minimum  for  power  law  interpolation  and  a  maximum 
for  logarithmic  interpolation.  The  difference  between  the  maximum  and  minimum 
Q-based  errors  in  Table  8  are  plotted  in  Figure  15.  Based  on  Figure  15  and 
Table  8,  we  henceforth  adopt  power  law  interpolation  (case  3,  Table  3)  as  the  pre¬ 
ferred  type  for  vertical  interpolations  of  Q.  Most  significantly.  Table  8  shows 
that  power  law  interpolation  of  Q  yields  processing  errors  that  are  very  competi¬ 
tive  with  the  smallest  errors  obtained  using  f  =  RH  or  f  =  Td> 


- F  =  LnO,  2-  LnP 

- F=0, Z=LnP 


Figure  15.  Vertical  Profiles  of  the  Global  Mean  BIAS  and  RMSE  for  the 
Synthesized  FGGE  III-A  Analysis  (17  January  1978,  00Z)  Obtained  Using 
Power  Law  Interpolation  of  Q  (Solid  Curve)  and  Logarithmic  Interpolation 
of  Q  (Dashed  Curve,  Repeated  From  Q-Based  Solid  Curve  in  Figure  9)  in 
the  NMC/AFGWC  Pre-  and  Postprocessor 


The  standard  logarithmic  interpolation  (F  =  f,  Z  =  LnP)  in  general  shows  the 
smallest  processing  errors  in  Table  8  for  the  RH-  and  T^-based  interpolations 
and  we  shall  retain  this  interpolation  type  henceforth  in  vertical  interpolations  of 
RH  or  T^.  However,  this  choice  was  less  clear  cut  than  that  for  Q  since  the  RH- 
based  and  T^-based  errors  in  Table  8  show  less  sensitivity  to  interpolation  type. 
To  explain  this  we  turn  to  Figure  16,  which  presents  the  T^  and  RH  profiles  cor¬ 
responding  to  those  for  Q  in  Figure  12(a).  Unlike  Q,  which  on  the  average  changes 
by  nearly  two  orders  of  magnitude  between  30  to  100  cb,  T^  and  RH  on  the  average 
change  by  less  than  25  and  50  percent,  respectively.  Hence  f  =  RH  or  f  =  T^  is 
less  responsive  than  f  =  Q  to  logarithmic  transformation. 
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Figure  16  shows  that  all  four  functional  choices  of  the  profile  are  quite 
similar  in  shape.  The  functions  F  =  T^,  Z  =  P  and  F  =  LnT^,  Z  =  P  are  both 
somewhat  concave.  Hence  in  Table  8,  both  yield  larger  bias  errors  (negative)  and 
somewhat  larger  RMSE  errors  than  the  functions  F  *  T^,  Z  =  LnP  and  F  =  LnTd, 

Z  =  LnP,  which  are  both  nearly  linear  and  which  both  yield  comparatively  small 
errors.  The  four  RH  profiles  in  Figure  16  are  also  similar  in  shape,  being  dis¬ 
tinctly  convex  at  70  cb  but  slightly  concave  at  85  and  40  cb.  Not  surprisingly  then 
in  Table  8,  the  standard  logarithmic  interpolation  for  RH  yields  a  fairly  large  pos¬ 
itive  bias  at  70  cb  and  slight  negative  biases  at  the  remaining  interior  levels.  The 
other  interpolation  types  for  RH  reduce  the  70  cb  positive  bias,  but  increase  the 
negative  bias  at  the  other  interior  levels.  Ultimately,  we  retain  logarithmic  inter¬ 
polation  for  RH,  because  when  averaged  over  the  six  mandatory  levels,  it  yields 
the  smallest  absolute  BIAS  and  RMSE. 

We  conclude  from  Table  8  that  logarithmic  interpolation  of  RH  or  T^  and  power 
law  interpolation  of  Q  all  yield  comparably  small  processing  errors.  We  will  adopt 
these  three  interpolation  schemes  as  our  nominal  methods  henceforth.  Table  9 
shows  the  sample  mean  processing  errors  and  their  sample  standard  deviation  over 
three  January  and  three  July  197  8  FGGE  III-A  analyses  for  the  nominal  interpola¬ 
tion  schemes.  In  Table  9,  the  results  shown  for  the  single  FGGE  III-A  analysis  of 
00Z,  17  January  1978  (used  exclusively  thus  far)  represent  well  the  results  of  the 
other  sample  dates.  The  differences  in  error  between  the  three  interpolation 
schemes  is  sufficiently  small  that  the  final  decision  to  pursue  one  particular  ap¬ 
proach  may  depend  on  secondary  factors  such  as  programming  convenience,  geo¬ 
graphical  distribution  of  the  error,  horizontal  truncation  error,  or  sensitivity  to 
spacing  of  a-layers.  We  shall  address  the  latter  three  factors  further. 

Regarding  sensitivity  to  o-structure,  we  suspected  that  the  close  proxi  iity  of 
the  <*3  layer  to  the  70-cb  mandatory  level  in  the  12 -layer  structure  of  Table  2  may 
be  suppressing  potentially  large  bias  errors  in  vertical  interpolations  using  RH. 

We  recall  that  the  hemispheric  mean  profile  of  RH  in  LnP  in  Figure  7  is  sharply 
convex  at  70  cb.  Therefore,  the  alternative  a- spacing  in  Table  10  was  tested. 

The  only  changes  in  Table  10  from  Table  2  occur  for  layers  Oy  and  o^,  which 
in  Table  10  are  further  separated  from  the  mandatory  levels  of  85,  70,  and  50  cb, 
respectively  (assuming  Pg  =  101.3  cb).  Figure  17  shows  the  sensitivity  of  the 
processing  error  to  the  new  o-spacing  for  the  three  nominal  interpolation  schemes. 

In  Figure  17,  the  logarithmic  interpolation  of  RH  indeed  yielded  the  largest 
increase  in  processing  error  for  the  new  o-structure,  with  Td  following  as  second. 
However,  the  Q-based  profile  was  notably  insensitive  to  the  change  in  o-spacing, 
especially  in  the  BIAS  error.  One  recognizes  from  Figure  11  that  the  more  linear 
is  a  given  function  F(Z),  such  as  F  =  LnQ,  Z  =  LnP  in  Figure  12,  the  less  sensi¬ 
tive  the  vertical  interpolation  error  will  be  to  the  location  of  the  o-layers.  We 
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Table  10.  The  12 -Layer 
CT-Structure  Used  in  the 
Pre-  and  Postprocessor 
Runs  for  Figure  17  (Com¬ 
pare  With  Table  2) 


expect  from  these  results  that  Q-based  power  law  interpolations  would  be  rather 
suitable  also  for  lower  vertical  resolutions  with  fewer  o-layers. 

All  (Sj,  A1>  processing  errors  shown  thus  far  have  been  vertical  profiles  of 
the  global  mean  BIAS  and  RMSE.  To  assess  the  geographical  distribution  of  pro¬ 
cessing  error,  we  show  in  Figure  19  the  latitude -height  contour  analyses  of  the 
zonal  mean  BIAS  and  RMSE  resulting  from  logarithmic  interpolation  of  RH,  Q, 
and  Tj  and  power  law  interpolation  of  Q.  The  zonal  averages  in  Figure  19  are 
based  on  the  input  winter  FGGE  III -A  analysis  of  00Z  on  17  January  1978  and  thus 
are  counterparts  to  the  global  mean  profiles  given  by  the  solid  curves  in  Figures 
9  and  15.  Also,  Figure  20  is  a  repeat  of  Figure  19,  but  for  the  input  summer 
FGGE  III-A  analysis  of  00Z  on  18  July  1978.  The  bottom  of  each  part  of  Figures 
19  and  20  shows  the  zonal  average  of  the  final  synthesized  analysis  for  the  respec¬ 
tive  interpolation. 

For  reference  purposes.  Figures  18(a)  and  (b)  show  the  zonal  average  of  the 
original  input  analyses  for  Figures  19  and  20,  respectively.  In  Figure  18,  the 
dry  tongues  of  the  descending  branches  of  the  Hadley  circulation  are  clearly  dis¬ 
cernible.  Moreover  the  Hadley  circulation  shows  a  reasonable  northward  shift 
from  January  to  July.  In  addition,  the  relative  humidity  decreases  with  height  at 
most  latitudes.  In  general,  then,  the  features  of  the  input  RH  analyses  appear 
reasonable.  We  caution,  however,  that  the  RH  in  the  upper-half  of  the  southern 
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Figure  17.  Vertical  Profiles  of  the  Global  Mean  BIAS  and  RMSE  for  the  Synthe¬ 
sized  FGGE  III-A  Analysis  (17  January  1978,  00Z)  Obtained  by  Executing  the  NMC / 
AFGWC  Pre-  and  Postprocessor  for  the  12 -Layer  cr-Structures  of  Table  2  (Solid 
Curve)  and  Table  10  (Dashed  Curve)  Using  Vertical  Interpolations  of  F  =  (RH,  LnQ, 
T<j)  in  Z  =  LnP 


hemisphere  troposphere  is  significantly  higher  than  that  in  the  corresponding  sea¬ 
son  of  the  northern  hemisphere.  Although  this  may  depict  true  inter -hemispheric 
variability,  it  also  is  very  likely  an  erroneous  artifact  of  the  forecast  model  in  the 
assimilation  system,  since  the  southern  hemisphere  is  persistently  data  sparse 
(especially  regarding  moisture  observations,  which  at  this  time  were  not  aug¬ 
mented  by  satellite  soundings  or  imagery  in  the  southern  hemisphere). 

In  the  processing  errors  in  Figure  19,  the  most  dramatic  result  is  the  lati¬ 
tudinal  pervasiveness  and  magnitude  of  the  large  errors  arising  from  logarithmic 
interpolation  of  Q  in  Figure  19(d).  Compared  to  the  other  plots  of  Figure  19,  this 
is  by  far  the  worst  choice  of  vertical  interpolation.  The  RH-based  logarithmic 
interpolation  in  Figure  19(a)  appears  to  give  the  most  satisfactory  processing  error 
over  all  latitudes  and  heights.  The  processing  errors  for  logarithmic  interpolation 
of  T^  in  Figure  19(c)  and  power  law  interpolation  of  Q  in  Figure  19(b)  are  also 


Figure  18.  Latitude -Height  Contour  Analysis  of  the  Zonal  Mean  of  the 
1978  FGGE  III -A  Analysis  of  RH  for  the  Winter  Case  (a)  and  the  Summer 
Case  (b)  Representing  the  Inputs  for  Figure  19  and  20,  Respectively. 

Coutour  interval  is  RH  =  10  percent 

small  and  competitive  with  those  in  Figure  19(a),  except  in  the  middle  troposphere 
at  high  latitudes  in  the  southern  hemisphere,  where  Figure  19(b)  and  19(c)  show 
larger  errors. 

The  above  discussion  of  Figure  19  also  holds  for  Figure  20.  We  further  note 
only  the  consistent  increase  in  RMSE  toward  both  poles  evident  in  parts  (a),  (b), 
and  (c)  of  both  Figures  19  and  20.  In  Figure  20,  this  poleward  increase  in  RMSE 
is  much  more  dominant  in  the  winter  hemisphere  than  in  Figure  19.  Further  in¬ 
spection  shows  this  poleward  increase  in  RMSE  is  not  the  result  of  a  closely 


39 


1 

81 

.  Z-LNP) 

JAN  1 

7. 1978 

OOZ 

S 

EQ 

8 

SYNTHESIZED  ANALYSIS 

(F«RH,  Z-LNP) 

JAN  17.1978  OOZ 

Figure  19.  Latitude -Height  Contour  Analysis  of  the  Zonal  Mean 
BIAS  (Top),  RMSE  (Middle),  and  Analyzed  RH  (Bottom)  for  the 
Synthesized  FGGE  III-A  Analysis  (17  January  1978,  00Z)  Obtained 
Using  Logarithmic  Interpolation  of  RH  (a),  Td  (c),  and  Q  (d),  and 
Power  Law  Interpolation  of  Q  (b).  Contour  interval  is  RH  =  2 
percent  (except  RH  =  10  percent  in  bottom  plot).  Dashed  contours 
denote  negative  values  (BIAS  only) 
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parallel  increase  in  the  BIAS  (see  RH-based  plots  in  particular).  However,  the 
next  section  on  spectral  truncation  error  will  provide  crucial  insight  into  the  cause 
of  the  maximum  of  RMSE  in  polar  regions. 

3.3  Horizontal  Truncation  Error 

In  Sections  3. 1  and  3.2  we  have  identified  and  eliminated  the  larger  sources  of 
bias  error  that  dominated  the  earlier  (S^,  A^)  error  statistics.  In  the  remainder 
of  Section  3,  we  shift  the  focus  primarily  to  RMSE  and  to  identifying  some  key 
sources  of  RMSE.  In  particular,  we  seek  to  determine  why  the  RMSE  consistently 
increases  toward  both  poles  in  Figures  19  and  20.  This  poleward  increase  of 
RMSE  occurs  at  all  levels  and  for  all  the  given  vertical  interpolation  approaches. 
Thus  we  speculate  that  its  cause  lies  with  horizontal  spectral  truncation  errors. 

Figure  21  from  top  to  bottom,  respectively,  shows  the  RMSE  plots  corres¬ 
ponding  to  parts  (a)  through  (c)  of  Figure  19,  but  wherein  all  the  spectral  trans¬ 
forms  in  the  pre-  and  postprocessor  [that  is,  steps  6  and  7  in  Table  5(a)  and  steps 
1  and  2  in  Table  5(b)]  have  been  skipped.  Clearly  the  poleward  increase  in  RMSE 
has  been  greatly  reduced  by  eliminating  the  spectral  transforms.  However,  much 
of  the  source  of  the  spectral  transform  error  lies  not  with  the  spectral  transform 
operation  itself,  but  rather  with  the  truncation  error  associated  with  the  particular 
moisture  field  being  transformed  --  namely  Q.  As  an  illustration,  in  Figure  22 
the  spectral  transforms  have  been  reinstated,  but  the  moisture  field  chosen  for 
spectral  transformation  in  each  plot  corresponds  to  that  used  for  vertical  inter¬ 
polation  (RH,  LnQ,  or  Td).  The  resemblance  between  Figure  22,  which  includes 
spectral  transforms,  and  Figure  21  is  dramatic.  Most  significantly,  neither 
Figure  22  nor  Figure  21  show  the  systematic  increase  of  RMSE  toward  the  poles 
that  was  present  in  the  RMSE  plots  of  Figure  19. 

The  moisture  variable  Q  is  substantially  more  prone  to  horizontal  spectral 
truncation  error  than  the  moisture  variables  of  RH,  LnQ,  and  Td  when  the  error 
is  measured  in  terms  of  equivalent  RH  error.  This  conclusion  is  highlighted  in 
Figure  23(a),  which  depicts  the  vertical  profiles  of  the  global  mean  RMSE  for  the 
three  spectral  transform  cases  just  discussed.  The  solid  curves  in  Figure  23(a) 
represent  the  standard  case,  repeated  from  the  bottom  of  Figure  17  and  corres¬ 
ponding  to  the  zonal  mean  RMSE  in  parts  (a)  through  (c)  of  Figure  19.  The  dotted 
and  dashed  curves  in  Figure  23(a)  correspond  to  Figures  21  and  22  respectively. 
The  results  in  Figure  23(a)  agree  with  the  results  of  Simmonds,  who  compared 
the  horizontal  spectral  truncation  errors  for  moisture  fields  of  RH,  Q  (or  mixing 

18.  Simmonds,  I.  (1975)  The  spectral  representation  of  moisture,  J.  Appl. 

Meteorol.  14:175-179. 
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Figure  21.  Latitude -Height  Contour  Analysis  of  the  Zonal  Mean 
RMSE  for  the  Synthesized  FGGE  III-A  Analysis  (17  January  1978, 
OOZ)  Obtained  by  Skipping  All  Spectral  Transforms  in  Pre-  and 
Postprocessing  While  Using  Logarithmic  Interpolation  of  RH  (Top), 
and  Td  (Bottom),  and  Power  Law  Interpolation  of  Q  (Middle). 
Contour  interval  is  RH  =  2  percent 
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Figure  23.  Vertical  Profiles  of  the  Global  Mean  BIAS  and  RMSE  for  the  Syn¬ 
thesized  FGGE  III -A  Analysis  (17  January  1978,  00Z)  Obtained  by  Executing 
the  NMC/AFGWC  Pre-  and  Postprocessor  Using  Vertical  Interpolation  of 
F  =  (RH,  LnQ,  or  Td)  in  LnP  With  the  Indicated  Choices  of  Horizontal  Spec¬ 
tral  Transforms  (a)  and  With  no  Transforms  for  Either  the  Original  or  Zero 
Terrain  (b) 


ratio),  and  dewpoint  depression  T-T^.  He  concluded  that  the  truncation  errors  of 
RH  and  T-T^  were  substantially  less  than  those  of  Q.  He  speculated,  as  actuaUy 
found  here,  that  T^  also  would  yield  less  truncation  error  than  Q. 

The  large  truncation  errors  found  here  near  the  poles  with  transforms  of  Q 
owes  to  the  fact  that  the  zonal  mean  value  of  Q  falls  off  rapidly  from  the  tropics  to 
the  poles  following  the  decrease  in  temperature.  Thus,  locally  small  spatial  de¬ 
viations  of  Q  in  the  tropics  may  lead  to  locally  large  deviations  in  Q  in  the  polar 
regions  following  spectral  reconstitution  of  the  grid  point  values  from  spectral 
coefficients.  Spectral  transform  theory  guarantees  only  minimization  of  the  global 
RMSE  between  input  field  and  spectrally  reconstructed  field.  This  global  minimi¬ 
zation  implicitly  forces  a  closer  fit  of  the  larger  Q  values  in  the  tropics  and  rela¬ 
tively  poorer  fit  of  the  smaller  Q  values  in  polar  regions.  In  the  cases  of  RH, 


LnQ,  and  T^,  the  decrease  in  zonal  mean  value  toward  the  poles  is  relatively 
much  less  (or  nonexistent  in  the  case  of  RH  —  see  Figure  18)  and  hence  the  stand¬ 
ard  deviation  of  these  fields  in  the  tropics  is  not  as  large  relative  to  the  respective 
zonal  mean  values  in  the  polar  regions  —  a  characteristic  that  yields  significantly 
smaller  spectral  truncation  error  in  polar  regions. 

The  standard  use  of  Q  as  the  moisture  variable  for  spectral  transformation  in 
the  pre-  and  postprocessor  here  is  of  course  dictated  by  the  use  of  Q  as  the  prog¬ 
nostic  moisture  variable  in  the  companion  GSM.  The  popularity  of  Q  as  the  prog¬ 
nostic  moisture  variable  in  many  global  forecast  models  stems  from  the  fact  that 
it  is  a  conservative  quantity  in  adiabatic  flow.  However,  the  cited  spectral  trunca¬ 
tion  errors  in  Q  in  spectral  forecast  models  make  it  difficult  to  maintain  the  phys¬ 
ical  constraint  that  Q  be  non-negative  and  casts  doubt  on  the  suitability  of  using  Q 
as  the  prognostic  moisture  variable  in  a  spectral  model.  Indeed,  based  on  Sim- 
monds'  work,  Daley  et  al*^  adopted  T-T^  as  the  prognostic  moisture  variable  in 
the  operational  GSM  of  the  Canadian  Meteorological  Center.  As  a  companion  argu¬ 
ment  against  using  Q  for  the  GSM  prognostic  moisture  variable,  we  propose  that 

3  3 

the  vertical  power  law  variation  of  Q  as  *  P  (that  is,  a  )  cited  in  Section  3.2 

20 

makes  it  particularly  vulnerable  to  vertical  truncation  errors.  Hansen  et  al 
briefly  discussed  the  vertical  differencing  of  Q  in  their  GCM  as  the  source  of  a 
systematic  positive  prognostic  moisture  bias.  As  an  apparent  solution  they  used  a 
harmonic  mean  in  place  of  the  arithmetic  mean  to  compute  values  of  Q  at  the  o- 
layer  interfaces  from  values  of  Q  at  the  o-layers.  In  a  planned  follow-on  study  to 
the  present  one,  we  hope  to  compare  the  suitability  of  various  moisture  variables 
including  Q,  LnQ,  T^,  and  T-T^  as  the  prognostic  moisture  variable  in  a  GSM. 

3.4  Terrain  Effects 

It  is  instructive  to  consider  the  consistent  increase  in  RMSE  at  100  cb  over 
that  at  85  cb  in  all  the  vertical  profiles  of  Figure  23(a).  This  increase  occurs  not 
only  for  all  the  given  types  of  vertical  interpolation,  but  also  irrespective  of  the 
presence  or  absence  of  spectral  transforms.  Tracing  the  cause  is  aided,  nonethe¬ 
less,  by  first  eliminating  the  spectral  transform  error  as  was  accomplished  in 
Figure  21.  Inspection  of  Figure  21  reveals  that  low-level  RMSE  is  a  maximum 
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at  middle  and  high  latitudes  in  the  N.  H.  and  at  high  latitudes  in  the  S.  H.  This 
pattern  suggests  a  relationship  with  latitudes  having  higher  zonal  mean  terrain 
heights. 

As  a  test  for  terrain  effects,  the  experiments  for  the  dotted  curves  in  Figure 
23(a)  were  repeated  wherein  a  globally  constant  terrain  height  of  1  m  was  imposed. 
The  resulting  global  mean  and  zonal  mean  RMSE  are  given  in  Figure  23(b)  and 
Figure  24  respectively.  These  figures,  by  their  absence  of  low-level  RMSE, 
verify  by  counterexample  that  the  presence  of  realistic  terrain  increases  low-level 
RMSE. 

The  terrain-induced  increase  in  100-cb  RMSE  stems  from  increased  low- 
level  extrapolation  errors  in  postprocessing.  We  recall  that  the  representative 
position  of  the  a-layers  given  in  Figure  7  assumes  a  surface  pressure  of 
Pg  =  101.3  cb.  Over  higher  terrain,  surface  pressures  will  be  substantially  less 
and  the  0-layers  will  be  displaced  substantially  upward  with  respect  to  the  stand¬ 
ard  pressure  surfaces.  Recalling  that  a 2  =  0.  862  from  Table  2,  we  find  that 
Pg  <  98.  6  cb  (that  is,  a  terrain  height  above  approximately  150  m)  implies 
P<02>  <  85  cb  and  thus  only  one  a-layer  will  exist  within  85  t  j  100  cb.  One  there¬ 
by  loses  the  exactness  of  the  low-level  linear  extrapolation  afforded  by  two  0- 
layers  occurring  within  85  to  100  cb,  as  depicted  in  the  right  of  Figure  8.  For 
surface  pressures  increasingly  less  than  98.  6  cb,  the  lower  extrapolation  error 
will  increase  as  the  extrapolation  interval  from  to  100  cb  increases.  Indeed 
the  larger  RMSE  at  100  cb  likely  is  dominated  by  the  higher  terrain  cases 

where  P  <  88.4  cb  (that  is,  terrain  height  above  about  1.  1  km),  when  the  lowest 
s 

a-layer  (a1  =  0.962)  is  above  the  85-cb  mandatory  level.  The  larger  100-cb 
RMSE  in  this  case,  however,  should  not  be  particularly  detrimental  to  the  con¬ 
tinuous  4-D  data  assimilation  process  as  the  P  -»  n  vertical  interpolation  in  the 
preprocessing  that  begins  each  GSM  first-guess  forecast  will  not  utilize  100-cb 
values  when  a^  is  above  85  cb. 

It  is  noteworthy  in  Figure  23(b)  that  although  the  presence  of  terrain  increases 
the  RMSE  at  100  and  85  cb,  it  slightly  decreases  the  RMSE  at  70,  50,  40,  and  30 
cb.  This  appears  consistent  with  the  fact  that  the  presence  of  terrain  raises  the 
mean  elevation  of  all  a-layers.  This  likely  displaces  layers  a y  a 4,  a5,  and  ag 
closer  on  average  to  their  nearest  standard  pressure  surfaces  (thus  decreasing 
the  interpolation  intervals)  than  that  shown  in  Figure  7.  The  latter  figure  assumed 
Pg  =  101.  3  cb,  which  is  representative  of  sea  level  (zero  terrain)  and  yields  a^, 
a4,  0g,  and  0g  layers  below  the  closest  mandatory  levels. 

An  earlier  portion  of  Section  3  discussed  how  verification  against  RH  analyses 
could  implicitly  favor  RH-based  vertical  interpolations.  Prior  to  the  no-trans- 
form  cases  of  Figure  23,  this  favoritism  in  RH-based  interpolations  was  largely 
negated  by  the  spectral  transform  steps,  which  required  conversion  from  RH  to  Q 
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Figure  25.  Repeat  of  the  Dot-Dashed  Curve  (a)  and  the  Solid  Curve  (b)  of 
Global  Mean  RMSE  in  the  RH-Based  Interpolations  of  Figure  23(b)  and  (a) 
Respectively,  Alongside  the  Parallel  Cases  (Dotted  Curves)  Where  Post- 
processed  RH  Values  are  the  Adjusted  Values  of  RH*  in  Eq.  (Bl) 


from  the  synthesized  temperature  field  in  a  manner  consistent  with  the  influence 
of  temperature  errors  in  the  Q-based  and  -based  interpolations.  The  develop¬ 
ment  and  discussion  of  Figures  25  and  26  are  not  crucial  to  the  conclusions 
reached  in  this  study  and  are  therefore  presented  in  Appendix  B.  The  reader  may 
choose  to  proceed  immediately  to  Section  4,  which  will  resume  discussion  of  the 
figures  with  Figure  27. 


4.  A  PSEUDO  4  D  DATA  ASSIMILATION 

A  major  concern  with  pre-  and  postprocessing  errors  lies  with  their  cumula¬ 
tive  impact  over  time  on  the  accuracy  of  an  intermittent  4-D  data  assimilation 
system.  As  discussed  in  the  introduction,  these  errors  are  among  many  sources 
of  error  that  can  contribute  substantially  to  undesirable  climatic  drift  and  error 
growth  in  the  assimilation  system.  As  a  means  of  measuring  the  error  impact  of 
moisture  pre-  and  postprocessing,  a  48-hr  "pseudo"  4-D  assimilation  was  exe¬ 
cuted.  That  is,  all  components  of  the  4-D  assimilation  in  Figure  1  were  run  for 
eight  cycles,  but  skipping  the  objective  analysis  stage  while  still  interrupting  the 
forecast  every  6  hrs  to  run  the  pre-  and  postprocessor.  The  resulting  pseudo 
assimilated  moisture  fields  were  verified  at  24 -hr  intervals  and  compared  with 
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Figure  26.  Latitude -Height  Contour  Analysis  of  the  Zonal  Mean  of  RH* 
in  Eq.  (Bl)  Corresponding  to  the  Global  Mean  of  RH*  in  Figure  25(a)  and 
Figure  25(b) 


the  verifications  of  a  parallel  standard  forecast  initialized  from  the  same  data 
(bottom.  Figure  1).  If  the  pre-  and  postprocessing  steps  in  the  pseudo  4-D  assim 
ilation  do  not  significantly  alter  the  moisture  fields  of  the  global  forecast  model, 
then  their  contribution  to  the  overall  error  growth  in  moisture  assimilation  can  be 
disregarded. 

To  measure  the  degree  to  which  this  is  the  case  for  this  pseudo  assimilation 
we  shall  extend  the  (F^,  A^)  and  (S^,  A^)  verifications  defined  at  the  start  of  Sec¬ 
tion  3.  The  verification  nomenclature  (F^,  A^)  and  (S^,  A^)  is  introduced  where 
the  former  denotes  the  verification  of  the  RH  fields  from  n  successive  executions 
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Figure  27.  Vertical  Profiles  of  the  Global  Mean  BIAS  and  RMSE  for  the  (Sj,  Aj), 
(Fb  ,  Ai),  and  (Fi,  A^)  Verifications  {Defined  in  Text)  Associated  With  the  48- 
hrTseudo  4-D  Assimilation  (n=8,  T=6  hr)  Starting  From  the  00Z  1978  FGGE  III- 
A  Analysis  of  17  January  [  19  January  for  (S?.  Aj)] .  Results  given  for  the  four 
vertical  interpolations  of  F  =  (RH,  LnQ,  Trf,  or  Q)  in  Z  =  LnP 


of  the  "update"  cycle  without  the  objective  analyses.  The  initial  FGGE  III-A  ob¬ 
jective  analysis  at  time  tQ  is  designated  Ao  while  the  verifying  analysis  at  time  t^ 
is  designated  A^.  Since  each  cycle  in  (F*JT,  A^)  is  T  =  6  hr  and  since  we  choose 
n  =  8,  the  overall  pseudo  assimilation  period  is  nT  =  48  hr.  The  quantity  (S^,  Aj) 
denotes  the  verification  of  the  RH  fields  from  n  successive  executions  of  only  the 
pre-  and  postprocessor  pair  applied  to  the  final  analysis  A^.  Thus,  (F"  ,  A^)  is 
a  dynamic,  temporally  evolving  forecast  from  the  initial  state  Aq,  whereas 
<S"  Aj)  is  again  a  static  transformation  or  synthesis  of  the  final  verifying  analy¬ 
sis  A..  Hopefully,  (S^  A.)  is  sufficiently  small  that  its  implied  effect  upon 
nT  ^  a  i  nT 

(Fj  ,  A^  is  minimal,  thereby  yielding  (F"  ,  Aj)  «  (F^,  A^),  where  F1  denotes 

a  standard  48-hr  forecast  execution. 

Figure  27  presents  results  of  a  pseudo  assimilation  starting  from  the  00Z 


FGGE  III-A  analysis  of  17  January  197  8.  The  global  mean  BIAS  and  RMSE  given 
for  the  (Sj.  A^),  (F”'*',  A^,  and  the  (F^,  A^)  verifications  are  valid  at  48  hr. 
Four  vertical  interpolation  alternatives  are  presented  from  left  to  right  in  Figure 
27  as  follows:  (1)  logarithmic  interpolation  of  RH,  (2)  power  law  interpolation  of 
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Q,  (3)  logarithmic  interpolation  of  T^,  and  (4)  logarithmic  interpolation  of  Q.  It 
should  be  emphasized  that  only  the  processing  steps  involve  the  alternative  mois¬ 
ture  variables,  while  all  forecast  steps  employ  Q  as  the  prognostic  moisture  var¬ 
iable.  The  forecast  model  used  in  these  experiments  is  the  AFGL  global  spectral 
model  (Brenner  et  al  )  with  the  same  R30  truncation  and  12 -layer  structure  as 
the  pre-  and  postprocessor.  Figure  28  shows  the  corresponding  temporal  trends 
of  the  global  mean  BIAS  and  RMSE  at  the  40-cb  level.  This  level  is  highlighted 
because  the  logarithmic  interpolations  of  Q  yield  the  largest  and  most  unacceptable 
errors  at  this  level.  Lastly,  Figure  29  and  Figure  30  show  the  zonal  mean  BIAS 
and  RMSE,  respectively,  corresponding  to  the  global  mean  plots  in  Figure  27  and 
in  the  same  order. 


Figure  28.  Temporal  Trend  at  P  =  40  cb  Level  for  the  Global  Mean  BIAS  and 
RMSE  of  the  (S^,  A^),  (F^Tf  Aj),  and  (F^,  Aj)  Verifications  in  Figure  27. 
The  (Sj.  Aj)  trend  is  shown  for  a  given  number  of  cycles,  n,  as  it  is  not 
strictly  defined  over  time  intervals 


21.  Brenner,  S. ,  Yang,  C. ,  and  Mitchell,  K.  (1984)  The  AFGL  Global  Spectral 
Model:  Expanded  Resolution  Baseline  Version,  AFGL-TR-84-0308. 
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Figure  29.  Latitude -Height  Contour  Analysis  of  the  Zonal  Mean  BIAS 
Valid  at  48  hr  (n  =  8,  T  =  6  hr)  and  Given  From  Top  to  Bottom  for  the 
(S^,  Ai),  (F9t,  Aj),  and  (Fi  ,  A j)  Verifications,  Respectively,  in  Fig 
ure  27.  Parts  (a)  through  (d),  respectively,  represent  the  interpola¬ 
tion  types  depicted  from  left  to  right  in  Figure  27 
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Looking  first  at  the  solid  BIAS  curves  in  Figure  27,  the  standard  48-hr  fore¬ 
cast  shows  a  modest  negative  bias  for  all  interpolation  types  except  at  100  cb. 
This  forecast  bias  is  significant  here  not  so  much  in  an  absolute  sense  but  rather 
in  defining  the  trend  of  the  forecast  when  devoid  of  processing  interruptions. 
Consistent  with  the  BIAS  plots  in  Figure  19,  the  dotted  BIAS  curves  in  Figure  27 
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show  that  the  synthesized  analyses  have  a  decidedly  positive  BIAS  for  all  interpo¬ 
lation  types  (drastically  so  for  logarithmic  interpolation  of  Q).  Not  surprisingly 
then,  the  dashed  BIAS  curves  for  the  pseudo  assimilations  (F^,  A^)  lie  between 

those  for  (S1?,  A.)  and  (F,,  A  ).  The  preferred  choice  of  vertical  interpolation  is 
li  li  nT 

that  which  yields  both  BIAS  and  RMSE  curves  for  (F^  ,  A^)  closest  to  those  for 
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the  associated  synthesized  analysis  is  greater  than  that  of  the  standard  forecast, 
a  result  not  even  closely  approached  by  the  other  interpolation  types. 

Recall  from  the  simple,  one -cycle  (S^  A^  error  statistics  of  Section  3.2 
that  Td -based  interpolations  yielded  competitively  small  processing  errors.  Con¬ 
trary  to  these  earlier  simple  tests,  the  present  T^-based  interpolations  over 
several  cycles  exhibit  an  appreciable  (though  not  overwhelming)  inferiority  to  the 
two  leftmost  interpolations  in  Figure  27.  Aside  from  the  logarithmic  interpola¬ 
tions  of  Q,  the  T^-based  interpolation  shows  the  largest  (S”,  A^)  BIAS  at  all  levels 
except  100  cb  and  consequently  the  greatest  spread  between  the  BIAS  curves  of  the 
(Fj,  A j)  and  (F"  ,  A^)  verifications.  Among  the  three  nominal  interpolations, 
we  conclude  that  the  BIAS  processing  error  of  T^-based  interpolations  is  more 
subject  to  continuous  accumulation  in  numerous  repetitions  of  pre-  and  postproces¬ 
sing.  This  conclusion  is  supported  by  Figure  28,  which  shows  that  of  the  three 
nominal  interpolations,  the  spread  between  the  (F“  ,  Aj)  and  (F^,  A^)  BIAS 
curves  at  40  cb  is  growing  most  rapidly  for  the  -based  interpolations. 

Thus,  logarithmic  interpolation  of  RH  and  power  law  interpolation  of  Q  re¬ 
main  as  the  most  likely  choices.  To  identify  a  preference  between  these  two  in¬ 
terpolations,  the  zonal  mean  BIAS  plots  in  parts  (a)  and  (b)  of  Figure  29  are  con¬ 
sidered.  At  70  cb,  both  interpolation  types  exhibit  a  disappointing  accumulation 
of  positive  BIAS  errors  in  their  (S”,  Aj)  verifications.  The  RH-based  interpola¬ 
tion,  though,  yields  a  somewhat  larger  BIAS  at  70  cb.  Further  evidence  of  this 
disadvantage  in  RH-based  interpolations  is  apparent  in  a  somewhat  poorer  preser¬ 
vation  in  part  (a)  over  part  (b)  of  the  (F^,  A^)  midtropospheric  BIAS  by  the 
(FjT,  ax)  BIAS. 

On  the  other  hand,  the  power  law  interpolation  of  Q  yields  a  much  larger 

accumulation  of  BIAS  errors  at  30  cb  in  the  top  of  part  (b)  than  in  the  top  of  part 

(a)  of  Figure  29.  The  influence  of  this  BIAS  error  is  clearly  present  in  the 
nT1 

(f"  ,  Aj)  BIAS  in  the  middle  of  part  (b).  Recalling  Figure  6,  the  underlying 
vertical  extrapolation  error  in  the  moisture  field  in  these  tests  is  virtually  zero 
at  30  cb  for  any  choice  of  moisture  variable.  Rather,  the  30  cb  BIAS  in  the  upper¬ 
most  plots  of  parts  (b)  and  (c)  of  Figure  29  were  traced  to  small  BIAS  errors  in 
the  (S”,  Aj)  temperature  verifications.  Thus,  rather  small  temperature  biases 
(that  is,  biases  in  saturation  Q  or  T^)  at  the  cold  temperatures  typical  of  30  cb 
can  cause  nontrivial  biases  in  derived  RH  values  computed  from  bias-free 
values  of  Q  or  T^.  This  effect  was  evidently  not  large  enough  to  make  RH  biases 
at  30  cb  a  dominant  feature  in  the  single  cycle  (S^,  A BIAS  plots  at  the  top  of 
Figure  19(b)  and  (c)  (although  the  initial  hints  of  the  problem  are  discernible). 

In  conclusion,  there  appears  to  be  no  clear  superiority  between  selecting 
logarithmic  interpolation  of  RH  or  power  law  interpolation  of  Q.  Either  of  these 
two  vertical  interpolations  is  suitable,  provided  in  all  cases  some  attention  is 


given  to  the  chosen  o-structure  (to  minimize  the  interpolation  interval  to  70  cb) 
and  to  the  magnitude  and  impact  of  temperature  biases  at  high  altitudes  (or  cold 
temperatures). 

5.  CONCLUSIONS 

This  study  examined  and  sought  methods  to  reduce  the  vertical  interpolation/ 
extrapolation  errors  of  moisture  in  the  pre-  and  postprocessor  of  an  intermittent 
4-D  data  assimilation  system,  such  as  that  used  at  AFGWC.  The  tests  of  vertical 
interpolation  utilized  the  three  moisture  variables  of  relative  humidity  (RH),  spe¬ 
cific  humidity  (Q),  and  dewpoint  (T^).  The  input  fields  for  the  pre-  and  post¬ 
processor  executions  were  chosen  from  a  sample  of  January  and  July  1978  FGGE 
III-A  objective  analyses. 

The  sample  mean  FGGE  III-A  vertical  profiles  of  RH,  Q,  and  T^  all  de¬ 
creased  with  height  at  all  moist  analyzed  standard  pressure  levels.  This  property 
dictated  that  the  constant  extrapolation  of  moisture  in  the  NMC/AFGWC  post¬ 
processor  be  replaced  with  linear  extrapolation.  This  is  particularly  true  when 
Q  or  Td  were  extrapolated.  When  the  prognostic  moisture  variable  is  carried 
only  at  model  cr-layers  in  the  troposphere,  one  can  further  reduce  upper  extrapo¬ 
lation  error  by  selectively  choosing  the  topmost  o-layer.  Since  30  cb  is  the  high¬ 
est  standard  pressure  surface  for  which  moisture  is  analyzed,  the  topmost  prog¬ 
nostic  layer  carrying  moisture  should  be  the  lowest  o-layer  consistently  above 
30  cb.  Lastly  regarding  extrapolation,  one  should  try  to  include  at  least  two 
o-layers  between  85  through  100  cb  to  further  minimize  lower  extrapolation  errors. 

Of  more  significance  are  the  results  obtained  for  vertical  interpolations  in¬ 
terior  to  the  modeled  moist  domain.  It  was  found  that  the  fairly  popular  use  of 
logarithmic  interpolation  (that  is,  linear  interpolation  in  LnP)  is  a  good  approach 
for  the  variables  RH  and  T^,  but  a  singularly  poor  approach  for  the  variable  Q. 
Logarithmic  interpolation  of  Q  leads  to  substantial  positive  moist  biases  in  the 
middle  and  upper  troposphere  because,  on  average,  Q  is  not  a  linear  function  in 
LnP.  Rather  it  varies  in  P  approximately  as  a  power  law  (that  is,  *  P  ).  Hence 
power  law  interpolation  of  Q  (that  is,  linear  interpolation  of  LnQ  in  LnP)  yields 
excellent  interpolation  accuracy.  In  arriving  at  the  above  result,  it  has  been 
demonstrated  that  it  is  good  practice  to  seek  an  elementary  functional  transforma¬ 
tion  F  =  F(f)  and  Z  =  Z(P)  of  the  moisture  variable  f  and  pressure  P  such  that  on 
the  average  the  function  F(Z)  is  closely  linear.  Straightforward  linear  interpola¬ 
tion  of  F  in  Z  can  thus  yield  interpolation  accuracy  greatly  surpassing  higher- 
order  interpolations  (for  example,  cubic  spline)  of  the  original  variable  f.  Also, 


vertical  interpolations  of  a  linear  function  F(Z)  are  less  sensitive  to  changes  in 
^-structure  or  resolution  and  submit  better  to  linear  extrapolation  when  neces¬ 
sary. 

The  vertical  interpolations  discussed  above  were  first  tested  in  a  static 
fashion  by  executing  just  the  NMC/AFGWC  pre-  and  postprocessor  pair.  The 
tests  were  then  extended  to  running  4  8 -hr  pseudo  4-D  intermittent  data  assimi¬ 
lations  wherein  the  objective  analysis  steps  in  Figure  1  (top)  were  eliminated 
while  still  interrupting  the  forecast  every  6  hr  to  execute  the  pre-  and  postproces¬ 
sor.  The  pseudo  assimilations  using  logarithmic  interpolation  of  Q  were  severely 
contaminated  by  the  positive  bias  discussed  above.  A  similar  exceptional  positive 
bias,  also  attributed  to  vertical  interpolations  of  Q,  was  previously  noted  by 
Lorenc  and  Tibaldi^  in  actual  ECMWF  4-D  data  assimilations  of  moisture. 

The  pseudo  assimilations  in  this  study  showed  that  logarithmic  interpolation 
of  RH  or  power  law  interpolation  of  Q  yield  comparably  good  results  in  minimiz¬ 
ing  the  accumulative  effects  of  vertical  interpolation  error  on  the  global  forecast 
model's  first-guess  fields.  Logarithmic  interpolation  of  T^,  though  performing 
well  in  the  static  tests,  showed  a  greater  tendency  toward  error  accumulation  in 
the  pseudo  assimilations.  Therefore,  it  is  not  to  be  recommended. 

It  was  noted  that  even  the  satisfactory  logarithmic  interpolation  of  RH  and 
power  law  interpolation  of  Q  both  yielded  nontrivial  positive  bias  tendencies  at 
70  cb  in  the  48-hr  pseudo  assimilations  [note  top  of  parts  (a)  and  (b)  in  Figure  29] . 
This  result  points  to  the  desirability  in  the  preprocessor,  even  when  using  fairly 
accurate  vertical  interpolations,  to  interpolate  the  analyzed  residuals.  This  in¬ 
creasingly  popular  approach  (for  example,  NMC,  ECMWF,  GLAS)  will  leave  the 
forecast  first-guess  fields  unmodified  the  update  cycle  in  data -void  regions. 
Nevertheless,  the  postprocessing  step  to  provide  forecast  first-guess  values  on 
pressure  surfaces  must  still  necessarily  interpolate  full -value  fields  from  a  to  P. 

There  are  other  alternative  4-D  data  assimilation  approaches  besides  use  of 
residuals  that  can  minimize  the  effects  of  vertical  interpolation  error.  For  ex¬ 
ample,  in  the  OI  objective  analysis  model  being  developed  at  AFGL  by  Norquist 
22  23 

et  al,  following  the  approach  of  Bergman,  the  analysis  is  performed  on  model 
n-surfaces,  eliminating  the  need  for  pre-  and  postprocessing.  It  does,  however, 
contain  a  one-way  interpolation  from  P  to  a  within  the  analysis  model  to  represent 
the  observations  on  cr-surfaces.  In  the  continuous  GFDL  4-D  assimilation  system 


22.  Norquist,  D. ,  Halberstam,  I.,  Tung,  S.  -L. ,  and  Johnson,  C.  (1984)  Two 

Methods  of  Global  Data  Assimilation,  AFGL-TR-84-0260. 

23.  Bergman,  K.  (1979)  Multivariate  analysis  of  temperature  and  winds  using 

optimum  interpolation,  Mon.  Wea.  Rev.  107:1423-1444. 
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described  by  Stern  and  Ploshay,  the  model  forecast  is  not  intermittently  inter¬ 
rupted.  Rather  OI -analyzed  observations  (full  values)  are  continuously  and  asyn- 
optically  inserted  into  the  model  forecast,  thereby  leaving  the  forecast  unaltered 
in  data -void  regions  as  well.  In  the  GFDL  assimilation,  however,  the  OI  analysis 
is  again  performed  on  pressure  surfaces  so  there  are  a -*  P  and  P  -*  o  interpola¬ 
tions  before  and  after  the  analysis. 

All  three  methods  cited  above  (residuals,  o-based  OI,  continuous  insertion), 
reduce  the  effects  of  vertical  interpolation  errors.  However,  they  still  include 
vertical  moisture  interpolations  at  key  points  in  the  assimilation.  Consequently, 
the  basic  sources  of  interpolation  error  demonstrated  in  this  study  will  still  im¬ 
pact  the  assimilations.  Therefore,  the  ease  with  which  the  moisture  interpola¬ 
tions  recommended  in  this  study  can  be  implemented  suggest  their  application  in 
all  data  assimilation  schemes. 

A  test  of  the  robustness  of  these  results  when  using  the  1979  FGGE  III -A 
analyses  of  NMC  is  presently  underway.  In  1979,  the  NMC  analyses  were  per¬ 
formed  using  optimal  interpolation  methods  rather  than  using  Hough  functions. 
Finally,  follow-on  investigations  are  underway  to  test  the  implications  of  the 
power  law  dependence  of  Q  in  the  truncation  errors  of  the  vertical  finite  differ¬ 
encing  of  Q-based  prognostic  moisture  equations  in  global  forecast  models. 
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Appendix  A 


Power  Law  Interpolation 


Suppose  given  quantity  f  satisfies  a  power  law  in  P  according  to 
f  =  aPb  ,  (Al) 


where  (a,b)  are  constants.  Given  known  values  of  f^  and  f 2  at  Pj  and  P^, 

Eq.  (Al)  can  be  applied  twice  to  obtain  the  two  simultaneous  equations  =  aPb 
and  f2  =  aPb.  Solving  these  equations  for  the  constants  (a,b)  yields 


b  = 


Lnf2  -  Lnf1 
LnP2  -  LnPj 


=  f  /Pb 

V  i 


Now  the  interpolated  value  of  f  at  P between  P^  and  P2  is  given  by 


f  =  aPD 


(A2) 


(A3) 


in  which  (a,b)  are  given  by  Eq.  (A2).  Finally  to  express  Eq.  (A3)  in  the  form 
of  Eq.  (1)  in  the  text,  the  logarithm  of  both  sides  of  Eq.  (A3)  is  taken  and  then 
Eq.  (A2)  is  substituted  for  (a),  which  after  some  algebra  gives 


Lnf  =  Lnfj  +  b  (LnP  -  LnPj) 


(A4) 


With  (b)  defined  in  Eq.  (A2),  the  above  result  is  equivalent  to  Eqs.  (1)  and  (2)  in 
the  text  for  the  case  of  F  =  Lnf  and  Z  =  LnP,  in  agreement  with  Table  3.  By  sim¬ 
ilar  analysis,  the  equivalence  of  Eqs.  (1)  and  (2)  to  the  remaining  interpolation 
types  in  Table  3  for  the  given  forms  of  F  and  Z  can  be  validated. 
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Appendix  B 


An  Alternative  Verification 


In  the  no-transform  cases  of  Figures  21,  23,  and  24,  which  eliminated  steps 
6  and  7  of  Table  5(a)  and  steps  1  and  2  of  Table  5(b),  the  pre-  and  postprocessing 
of  the  input  moisture  field  (RH,  LnQ,  or  T^)  proceeded  fully  independent  of  the 
temperature  field  until  the  final  conversion  to  RH  before  verification  [  step  4  in 
Table  5(b)].  In  RH-based  interpolations,  even  this  last  step  was  unnecessary. 
Thus,  in  the  no-transform  RH-based  interpolations,  the  processed  moisture  field 
was  completely  uncontaminated  by  temperature  processing  errors.  For  example, 
the  accuracy  of  the  upper  and  lower  extrapolation  of  the  moisture  field  in  the  no¬ 
transform,  zero-terrain  RH-based  interpolations  is  particularly  good  at  100  and 
30  cb  (virtually  zero  RMSE)  at  the  left  of  Figure  23(b)  and  top  of  Figure  24.  This 
follows  because  as  the  right  of  Figures  6  and  8  showed,  over  zero  terrain  the 
present  extrapolation  procedures  exactly  reproduce  the  input  moisture  value  f  at 
100  and  30  cb.  In  the  no-transform,  zero -terrain,  LnQ-  and  T^-based  interpola¬ 
tions  of  Figures  23(b)  and  24,  the  lower  and  upper  extrapolation  of  the  moisture 
variable  (LnQ  or  T^)  is  still  exact  but  now  obscured  by  the  verification  step, 
which  as  stated  earlier  requires  postprocessed  temperature  fields  to  convert 
postprocessed  LnQ  or  T^  to  RH.  In  all  three  no-terrain  profiles  of  Figure  23(b), 
the  corresponding  global  mean  RMSE  of  temperature  was  0.  17,  0. 19,  0.21,  0.34, 
0.28,  and  0.33  K  at  100,  85,  70,  50,  40,  and  30  cb  respectively. 

To  impose  a  representative  effect  from  processed  temperature  error  on  the 
verification  of  the  no -transform ,  zero-terrain,  RH-based  interpolation  of  Figure 
23(b),  we  modified  the  postprocessed  RH  field  with  the  adjustment  given  by 
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where  (~)  denotes  a  postprocessed  value  and  T  denotes  the  original  input  temper¬ 
ature  value.  To  motivate  the  chosen  adjustment  in  Eq.  (Bl),  we  note  that  in  the 
LnQ-  or  T  -based  interpolations,  step  4  of  Table  5(b)  involves  first  obtaining  e 
directly  from  Q  or  T  ,  (as  yet  unaffected  in  no-transform  cases  by  processing 
errors  in  temperature)  and  finally  dividing  by  e  (T)  to  obtain  RH  according  to 

Eq.  (6).  Likewise,  RH  X  e  (T)  in  Eq.  (Bl)  represents  an  expression  for  e  that  is 

s 

unaffected  by  temperature  processing  errors,  which  are  lastly  imposed  by  the 
division  by  e  (T). 

°  <v 

The  comparison  of  the  RMSE  verifications  of  RH  and  RH*  in  the  no-transform, 
zero-terrain,  RH-based  interpolation  is  given  in  Figure  25(a).  The  corresponding 
zonal  mean  cross-section  of  RMSE  for  RH*  is  given  in  Figure  26(a)  (compare  with 
top  plot  in  Figure  24).  As  expected,  the  RMSE  at  100  and  30  cb  in  these  figures  is 
noticeably  larger  for  RH*  than  for  RH.  Clearly,  the  RMSE  of  RH*  in  Figures 
25(a)  and  26(a)  is  a  fairer  measure  than  that  of  RH  for  comparison  with  the  zero- 
transform,  no-terrain,  LnQ-  and  Td-based  RMSE  results  of  Figures  23(b)  and  24 
respectively.  The  close  resemblance  between  Figure  26(a)  and  the  lower  two 
plots  of  Figure  24  is  particularly  noteworthy  (especially  in  the  N.  H. )  and  serves 
to  highlight  once  again  the  nearly  equal  performance  among  the  choices  of  RH, 

LnQ,  or  Td  as  the  moisture  F  variable  in  vertical  interpolations  in  Z  =  LnP. 

With  the  above  as  background,  we  next  examine  the  impact  of  the  favoritism 
toward  RH-based  interpolations  in  the  verification  of  the  standard  case  that  in- 
eludes  full  transforms  in  Q  and  realistic  terrain.  Figure  25(b)  compares  the  RH- 
derived  and  RH*-derived  profile  of  global  mean  RMSE  and  Figure  26(b)  shows  the 
ffH*-derived  cross  section  of  zonal  mean  RMSE  [compare  with  -derived  RMSE 
plot  of  Figure  19(a)].  As  anticipated,  the  optimism  toward  RH-based  interpola¬ 
tions  in  verifying  the  standard  cases  with  full  transforms  in  Q  and  realistic  ter¬ 
rain  is  negligible  at  all  levels,  except  at  100  cb.  Even  for  the  latter,  the  cross 
section  of  zonal  mean  RMSE  in  Figure  26(b)  vs  that  in  Figure  19(a)  shows  the  bulk 
of  the  increased  RMSE  for  adjusted  RH*  to  be  present  at  the  highest  latitudes. 
Lastly  and  most  significantly,  comparison  of  Figure  26(b)  with  the  RMSE  plots  of 
Figures  19(b)  through  (d)  upholds  the  conclusion  of  Section  3.2  that  the  (S ^ ,  Aj) 
error  characteristics  for  vertical  interpolations  in  LnP  using  RH,  LnQ,  or  Td 
are  very  competitive  with  one  another  (and  all  vastly  superior  to  vertical  inter¬ 
polations  of  F  =  Q  in  Z  =  LnP). 


LIST  OF  ACRONYMS 


AFGL  Air  Force  Geophysics  Laboratory 

AFGWC  Air  Force  Global  Weather  Central 

AWAPS  Advanced  Weather  Analysis  and  Prediction  System 

ECMWF  European  Center  for  Medium  Range  Weather  Forecasts 

FGGE  First  GARP  Global  Experiment 

FNOC  Fleet  Numerical  Oceanography  Center 

GARP  Global  Atmospheric  Research  Program 

GFDL  Geophysical  Fluid  Dynamics  Laboratory 

GLAS  Goddard  Laboratory  for  Atmospheric  Sciences 

GSM  global  spectral  model 

HIRAS  High  Resolution  Analysis  System 

NMC  National  Meteorological  Center 

NMI  normal  mode  initialization 

NWP  numerical  weather  prediction 

OI  optimal  interpolation 
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